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SUMMARY 


Studies  of  nitrogen  transport  In  materials  has  been  limited 

by  the  lack  of  a long-lived  nitrogen  radioisotope.  Of  the 

seven  known  nitrogen  isotopes  only  two  are  found  in  nature, 

and  both  are  stable.  The  short  half-lives  of  the  unstable 

isotopes  have  prevented  their  utilization  in  the  standard  radio- 

tracer  experiments.  This  study  utilized  the  direct  observation 
15  12 

of  the  N(p,c(q)  C nuclear  reaction  to  obtain  nitrogen  self- 

diff usivities  in  both  the  grains  and  grain  boundaries  of  NOR- 
TM 

ALIDE  NC-132  hot  pressed  silicon  nitride.  The  values  of  the 
self-diffusion  coefficients  measured  were: 


2.02  x 10 


-11 


.+4.90  x 10 
-2.01  x 10 


-9 

-ID 


exp 


(3.17  - 0.74)  x 105 
RT 


in  the  grains,  and 

DGBfc>  “ *-°9  <-42:039  X ^^exp  - LI -09^  1.72)  »_lol 

RT 

in  the  grain  boundaries. 

A discussion  of  nuclear  microanaly tical  techniques  and  their 

utlization  in  the  field  of  materials,  and  the  calculation  and 

15  12 

measurement  of  the  N(p,aQ)  C differential  cross  section  at  a 
laboratory  scattering  angle  of  165°  is  also  included. 
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SECTION  I 
INTRODUCTION 

Silicon  nitride,  Si^,  one  °f  the  most  interesting  ceramic 
materials  to  be  developed  in  recent  years,  has  attracted  much  atten- 
tion due  to  its  possible  utilization  as  a high-temperature  material 
for  gas  turbines,  bearings,  and  other  advanced  structural  applica- 
tions. Silicon  nitride  has  good  strength  in  its  hot  pressed  form, 
a low  coefficient  of  thermal  expansion,  and  a high  thermal  conduc- 
tivity which  combine  to  give  it  a virtually  unsurpassed  thermal  shock 
resistance.  While  it  does  oxidize  readily,  its  oxidation  resistance 
is  relatively  good  due  to  the  adherent  silica-based  layer  formed. 

To  date  most  of  the  information  (1)  available  about  silicon 
nitride  deals  with  the  production  of  a material  with  optimum  mechani- 
cal properties  or  with  the  microstructure.  In  order  to  better  un- 
derstand the  processing  kinetics  of  Si3N4,  detailed  measurements  of 
nitrogen  transport  are  valuable,  but  incomplete  to  date. 

Nitrogen  transport  studies  in  nitrides  have  been  limited  by  the 
lack  of  isotopic  tracers.  Most  measurements  (see  Chapter  III)  have 
been  inferred  from  reaction  kinetic  data,  while  a few  studies  have 
followed  the  decrement  in  isotopic  concentration  of  a gas  in  contact 
with  a solid  sample.  This  study,  to  the  author's  knowledge,  is  the 
second  application  of  a nuclear  technique  to  the  study  of  nitrogen 
transport,  and  the  first  application  of  the  direct  observation  of  a 
nuclear  reaction. 


Because  the  author  has  a considerable  interest  in  nuclear 
microanalysis  and  nuclear  physics.  Appendices  are  devoted  to 
these  topics. 


SECTION  II 

DIFFUSION  PHENOMENOLOGY 


2.1  General  Laws 

For  the  case  considered  here,  diffusion  can  be  regarded  as  a 
concentration-leveling  process.  Tracer  diffusion  involves  the  ex- 
change of  a marked  and  unmarked  species  and  hence  at  least  two  dif- 
fusion equations  can  be  considered.  However,  due  to  coupling  of  the 
equations  only  one  species  needs  to  be  looked  at  explicity. 

The  mathematical  basis  for  the  diffusion  process  is  found  in 
Fick's  laws.  Fick's  First  Law  relates  the  flux  of  diffusing  spe- 
cies,  J to  the  concentration  gradient  of  the  diffusing  species,  vC, 
(where  C is  in  atom  fraction)  by  the  following  equation: 

S = -pDvC  (1) 
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where  D is  the  diffusion  coefficient  and  p is  the  density.  Often 
one  is  concerned  with  diffusion  in  one  direction  only,  and  Equa- 
tion 1 reduces  to: 


J = -pD(lC) 

X 3X 


(2) 


Equation  1 and  2 agree  with  the  empirical  fact  that  when  the  con- 
centration gradient  goes  to  zero,  the  driving  force  goes  to  zero, 
and  there  are  no  external  forces  acting  on  the  sample,  the  net  flux 
of  diffusing  species  goes  to  zero. 

Fick's  Second  Law  relates  the  time  rate  of  change  in  concen- 
tration with  the  flux  by  the  following: 


/jjC\  = - vj  = v(DvC) 
v3t'x,y,z  t 


(3) 


and  reduction  to  one  dimension  yields: 


<§)„  - |-  (D 
3t  x 3X  3x't 


(4) 


If  the  diffusion  coefficient  is  concentration  independent.  Equation  4 
can  be  rewritten  as: 


(|f)x  s 0<^>t 


3x 
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(5) 


The  solution  of  Equation  5 will  depend  on  the  initial  and 
boundary  conditions  imposed  by  the  particular  experiment.  For 
the  case  where  a semi-infinite  solid  is  exposed  to  a gas  of  con- 
stant composition  and  there  is  no  phase  boundary  reaction  at  the 
surface  the  following  initial  and  boundary  conditions  apply: 

C(x.O)  = CQ  (6) 

which  states  that  the  initial  concentration  of  diffusing  species 
is  constant  throughout  the  sample, 

C(-,t)  = Cq  (7) 

states  that  the  sample  is  a semi-infinite  medium,  and 

C(0,t)  « Ce  (8) 

The  concentration  of  diffusing  species  at  the  sample  surface  is 
constant  with  time.  The  solution  to  the  differential  equation 
is  given  by: 


C(x.t)  * C + (C  - CJ  erfc  ( x ) 

o e o ^ (9) 

2(Dt) 
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where  erfc  ( is  the  complementary  error  function  of  ( — - — r) 

2(Dt)  2(Dt) 

and  is  given  by: 


2(Dt) 


erfc  (-5—  ) = 1 - L- 

2(Dt) 


-w 

e dw 


The  temperature  dependence  of  the  diffusion  coefficient  is 
usually  given  by: 


D = D exp  - JB* 
o RT 


where  Dq  is  a frequency  factor,  Q*  is  the  activation  energy  for 
diffusion,  R is  the  gas  constant,  and  T is  the  temperature. 


2.2  Diffusion  with  a Finite  Rate  of  Surface  Exchange 


The  situation  often  arises  where  it  is  necessary  to  relax  the 
boundary  condition  given  by  Equation  8,  that  the  concentration  of 
diffusing  species  at  the  sample  surface  is  a constant  with  time,  to 
allow  for  a finite  rate  of  exchange  at  the  surface.  The  presence  of 


r 


a surface-exchange  coefficient,  K$,  which  is  a first-order  reac- 
tion-rate constant  will  effect  both  the  shape  and  magnitude  of 
C(x,t).  When  a surface-exchange  coefficient  is  required,  the 
boundary  condition  given  by  Equation  8 is  written  as: 


Ks(Ce  - C(0,t)  = - 


The  solution  to  Fick's  Second  Law,  Equation  5,  for  this  new  set 
of  initial  and  boundary  conditions  has  been  given  by  Robin  (2)  as: 


C(x,t)  = CQ  - (Ce  - C0)  {exp( 


K (x  + Kt 
-s fT-5-)  erfc 


x + 2K  t 

r s. 


2(Dt) 


- erfc  (— x u)  } 

2(Dtr 


As  can  be  seen  if  K$  -*• » Equation  13  reduces  to  Equation  9. 


2.3  Diffusion  in  Bicrystals 


In  polycrystalline  materials  high  diffusivity  paths,  grain 
boundaries,  exist  within  the  sample.  One  should  therefore  expect 
the  concentration  profiles  obtained  to  be  due  to  diffusion  along 
grain  boundaries  and  diffusion  into  the  grains  from  both  the  sur- 
face and  the  grain  boundaries. 
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Solutions  to  the  grain  boundary  diffusion  problem  have  been 
presented  by  Fisher  (3),  Whipple  (4),  and  Suzuka  (5)  and  analyzed 
by  LeClaire  (6).  The  configuration  used  in  the  mathematical  treat- 
ment of  grain-boundary  diffusion  is  shown  schematically  in  Figure  1. 
The  grain  boundary  is  treated  as  a slab  of  uniform  composition  and 
thickness  6.  It  is  assumed  that  Dgg>>  Dy.  Fick's  Laws  can  now  be 
solved  in  the  grain  boundary  and  within  the  grains  by  applying  the 
appropriate  initial  and  boundary  conditions.  The  case  of  a constant 
surface  concentration  has  been  treated  by  both  Fisher  (3)  and  Whipple 
(4)  and  will  be  examined  here.  The  initial  and  boundary  conditions 
that  apply  are  given  by: 


C(x,y,0)  - CQ  = 0 

(14) 

c(-,y,t)  - c0  = 0 

(15) 

C(0,y,t)  = Ce 

(16) 

Continuity  of  both  concentration  and  flux  must  also  be  maintained 
at  the  interface  between  the  grain  boundary  and  grain.  These  continuity 


I 

i 
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conditions  are  given  by: 


Cy(x,  t ~ 6,t)  - CGB(x,  ± j 6,t) 
3Cv  3CGB 

Dv(aP_)J56  * dgb(  sy 


(17) 


(18) 


Within  the  grain  boundary  C will  be  an  even  function  of  y and 

2 

can  therefore  be  approximated  by  a power  series  in  y with  the  con 

stants  being  functions  of  x and  t. 

From  the  boundary  conditions  and  Fick's  Second  Law  the  follow- 

2 

ing  can  be  obtained,  ignoring  terms  of  order  6 : 


2D 


/1C) 


r32C 

+ °GB(3  2^6 
x 


3C 

3t 


(19) 


which  must  be  solved  at  the  boundary.  The  first  term  represents 
diffusion  out  of  the  grain  boundary  and  the  second  term  represents 
diffusion  along  the  grain  boundary.  The  equation  to  be  solved  within 
the  grains  is  given  by: 


d A)  * 

v ax2  ay2  c 


(1£) 
at  X,y 


(20) 


Following  the  derivation  of  Whipple  (4)  it  is  convenient  to  con- 


vert to  dimensionless  variables,  defined  by: 


Fisher  (3)  has  obtained  the  following  approximate  solution  for 


C = C exp(-u4  y t *S)erfc(4  c)  (24) 

e / 

by  assuming  that  y is  sufficiently  large  and  that  only  diffusion  from 
the  grain  boundary  contributes  to  C. 


Whipple  (4)  has  provided  the  exact  solution  to  Equation  20  pro- 
°GB 

viding  y > — — , which  is  usually  the  case.  The  exact  solution  is: 


where  w is  a dummy  variable  of  integration.  The  first  term  in 
Equation  25  is  the  direct  volume  diffusion  contribution  to  C and  is 
the  one -dimensional  solution  to  the  diffusion  equation.  The  second 
term  gives  the  contribution  to  C from  diffusion  out  of  the  grain 
boundaries. 

If  t >1,  the  second  term  in  Equation  25  can  be  rewritten  as  (4): 


Equation  26  is  known  as  Whipple's  asymptotic  solution.  If  the  quantity 

-h 

YT  >>  1 then  the  first  term  in  Equation  25,  that  due  to  direct  volume 
diffusion,  can  be  neglected  and  C will  be  given  by  Equation  26. 

From  Equations  21  and  23  the  following  can  be  obtained: 


i 

I 

» 

I; 


dgbs  ■ ) "2  (~r>h  >2 

GB  3X  C 3(Yt  *> 


The  bar  indicates  average  concentration  with  respect  to  y;  and  it  is 

-V 

understood  that  differentiation  with  respect  to  yr  is  with  t and  Dy 
held  constant. 

The  last  term  of  Equation  27,  (^-lnC  , ) is  a constant  as  can  be 

9(yr) 

seen  by  integration  of  the  Fisher  solution,  Equation  24,  with  respect 
to  c(6): 


-u 

Here  the  last  tern  Is  simply  n and  Is  Independent  of  C and  x. 


Similarly,  It  Is  possible  to  obtain  an  expression  for 
3 InC  2 

( j )c  from  Whipple's  asymptotic  solution.  Equation  26: 

3 (yt  S) 


,31nC  .2 

3 (yt  *) 


- [-  f (yt_ 1 - 0.6307  (yt-*5)173 

+ 0.264  ] 

T 1 


(29) 


1)  It  should  be  noted  that  Equations  28  and  29  apply  only  to 
regions  where  direct  volume  diffusion  can  be  neglected.  2)  The 

— Jj 

asymptotic  solution  is  a good  approximation  if  the  value  of  yt  >9. 
3)  The  asymptotic  solution  is  not  a constant.  However  for  most  ex- 

-\c 

perimencal  data  the  midpoint  of  the  range  of  yt  can  be  used  due  to 
the  low  curvature.  4)  Equation  27  will  give  values  of  D..d5  approxi- 
mately  three  times  too  small  if  the  Fisher  solution  is  used  and 
approximately  10%  too  large  if  Whipple's  asymptotic  solution  is  used 

-jj 

when  yt  falls  between  4 and  8(6). 

3 InC 

Levine  and  MacCallum  (7)  noted  that  the  quantity  r A/r 

3(yt'V/5 

_L 

was  essentially  independent  of  (yt  ) over  the  approximate  range 

-V 

2 < YT  ^ < 10.  The  quantity  DGg^  can  then  be  calculated  from: 


(30) 

which  is  analogous  to  Equation  27.  For  large  values  of  t(t  > 10)  the 

value  of  is  a constant  equal  to  0.78. 

3(yt 

Borisov  and  Lyubov  (8)  report  this  ratio  to  have  a value  of  0.99 

~h 

for  values  of  (yx  ) < 1. 

Lundy  and  Federer  (9)  presented  a graphical  method  for  obtaining 
the  grain-boundary  and  volume  diffusivities  given  the  concentration 
versus -depth  profile.  The  logarithm  of  the  concentration  is  first 
plotted  as  a function  of  x(x^"*  if  the  model  of  Levine  and  MacCallum 
is  used).  A straight  line  is  then  fit  to  the  points  corresponding 
to  the  deepest  penetrations.  The  slope  of  this  line  can  be  used  in 
Equation  27  (or  Equation  30)  to  determine  the  value  of  D 6 . The 
concentration  distribution  corresponding  to  the  linear  fit  is  then 
subtracted  from  the  measured  concentration  profile.  The  resulting 
profile  is  then  due  to  volume  diffusion  only  and  can  be  fit  with  the 
appropriate  solution  to  Fick's  Second  Law. 

2.4  Extension  to  the  Polycrystalline  Case 

The  equations  derived  in  the  previous  section  are  for  diffusion 
in  a bicrystal  however  we  wish  to  apply  them  to  polycrystalline 
materials.  It  Is  therefore  of  Importance  to  examine  what  meaning 


31nC 


GB  3x6/5  C 3(Yt 


>5/3 


-y/s' 


the  calculated  dif fusluities  will  have  when  these  equations  are  ap- 
plied to  polycrystals.  If  it  is  assumed  that  the  grain  size  is  a, 
then  this  value  can  be  related  to  the  characteristic  volume  diffusion 

j, 

distance,  2(Dvt)  , and  the  characteristic  grain  boundary  diffusion 

distance,  2(D_Bt)  . 

GB 

v I u 

The  first  possible  case  is  when  2(D  t)^  < — a and  2 (D  t)^  < a. 

v 2 GB 

In  this  case  the  value  of  D and  D will  be  correct  if  all  the  grain 

boundaries  are  perpendicular  to  the  surcase.  However,  as  is  the 

case  for  real  materials  the  grain  boundaries  intersect  the  surface  at 

all  angles.  Use  of  Equation  27  or  30  will  lead  to  an  underestimation 

2 

of  Dgb  by  factor  of  — . If  the  plane  of  the  grain  boundary  makes  an 

$ with  the  surface  of  the  sample,  see  Figure  2a,  the  measured 

penetration  distance  in  the  x direction  will  be  given  by  b sin  $. 

As  can  be  seen  if  <t>  is  90°  x will  equal  b and  if  $ is  0°  or  180° 

x will  be  0,  that  is,  no  penetration  will  be  seen.  According  to 

Equations  27  and  30  the  calculated  value  of  D will  be  directly  pro- 

GB 

2 

portional  to  (3x)  . To  determine  the  amount  by  which  D is  under- 

GB 

estimated  for  a uniform  distribution  of  grain  boundary  angles  the 

2 

average  value  of  (b  sin  $)  must  be  found  over  a hemisphere: 


«b  sin  »)2>  - L&-0 


(31) 


Because  the  process  of  subtracting  the  values  associated  with  either 


(—  ) or  (^575)  from  the  experimental  data  removes  the  contri 

dx  ' 


Figure  2a:  Configuration  of  a tilted  grain  boundary. 


Figure  2b:  First  possible  configuration  for  intersecting  grain 
boundaries . 


Figure  2c:  Second  possible  configuration  for  intersecting  grain 
boundaries . 
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tD  the  cd rce rtratio n versus  depth  profile  from  both  diffusion 
along  grain  boundaries  and  diffusion  into  the  grains  from  the  grain 
boundaries  the  resulting  profile  should  be  due  to  diffusion  from  the 
surface  only.  Therefore  the  value  of  calculated  will  not  be  signi- 
ficantly different  from  the  actual  value. 

The  second  possiblity  involves  tilted  intersecting  grain  boun- 
daries, i.e.  2(D  t)Jj>  a while  2(D  t)  i a.  The  fact  that  the 
Gl>  V * 

boundaries  are  tilted  will  lead  to  an  underestimation  of  D^,  because 
the  actual  penetration  distance  along  the  grain  boundaries  will  ap- 
pear to  be  less  than  if  the  boundaries  were  perpendicular  to  the  sur- 
face. To  determine  the  amount  by  which  DGB  is  underestimated  It  is 
necessary  to  examine  the  two  types  of  intersecting  boundaries  en- 
countered. The  first  type  has  one  boundary  branching  to  form  two  new 
boundaries,  see  Figure  2b.  In  order  to  maintain  a diffusing  species 
flux  balance  at  the  intersection  the  flux  along  grain  boundaries  II 
and  III  will  be  less  than  that  along  grain  boundary  I.  The  second 
type  of  intersection  is  for  two  boundaries  to  come  together  to  form 
a third  as  in  Figure  2c.  In  this  case  to  maintain  a flux  balance  at 
the  intersection  the  flux  along  boundary  III  will  be  greater  than 
the  flux  along  either  boundary  I or  II.  Statistically,  the  number 
of  each  type  of  intersecting  boundaries  should  be  the  same,  thus  the 
net  effect  of  the  intersections  on  the  calculated  value  of  should 
not  be  significant.  Thus,  D.„  will  be  underestimated  by  approximately 

uD 

2 

. For  the  reasons  given  in  the  first  case  Dv  will  not  be  signi- 
ficantly affected. 


■i  i - 
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The  final  possibility  also  involves  tilted  intersecting  grain 
W u i 

boundaries,  2(DGfit)  > a,  but  (Dvt)  > y a.  For  reasons  identical 

2 

to  those  above  DGfi  will  be  underestimated  by  approximately  • The 
diffusion  profiles  from  opposite  sides  of  the  grains  will  now  over- 
lap, thus,  decreasing  the  concentration  gradient.  This  will  tend 

to  overestimate  D . 

v 


Section  III 

PREVIOUS  STUDIES  OF  NITROGEN  TRANSPORT 

3.1  General 

The  possibility  of  studying  nitrogen  transport  in  materials  has 
been  limited  by  the  availability  of  isotopic  tracers.  Of  the  seven 
known  nitrogen  isotopes,  only  two  are  found  in  nature  and  both  are 
stable,  see  Table  1.  The  short  half-lives  of  the  unstable  Isotopes 
prevents  their  utilization  in  usual  radiotracer  experiments.  Nitro- 
gen diffusion  measurements  have  therefore  been  confined  to  gaseous 
exchange,  studied  by  following  the  change  in  isotopic  concentration 
of  the  gas,  activation  analysis,  or  inferred  from  reaction  data. 

Smeltzer  and  Desmaison  (10)  in  an  appraisal  of  the  literature 
on  the  measured  diffusion  coefficients  of  nitrogen  in  metal  nitrides 
conclude  that  the  studies  remain  at  an  infancy  state.  Table  2 (10) 
gives  the  available  diffusion  coefficients  for  some  binary  nitrides. 

The  scatter  in  both  the  measured  frequency  factor,  Dq,  and  activation 

* 

energy,  Q , can  be  attributed  to  the  lack  of  suitable  single  cry- 
stals and  the  lack  of  suitable  techniques  to  measure  nitrogen  trans- 
port. 

Most  studies  have  been  carried  out  on  polycrystalline  materials 
which  often  have  unknown  structures  and  impurity  levels.  This  has 
led  to  difficulties  in  separating  the  volume  diffusion  coefficient 
from  the  contributions  due  to  grain  boundaries,  dislocations,  and 
Impurity  effects. 


i 
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TABLE  1 

KNOW  NITROGEN  ISOTOPES 


Isotope 

X Abundance 

Decay  Mode 

Half-Life 

l> 

+ 

6 

0.011  Sec 

+ 

B 

599.4  Sec 

> 

99.63 

0.37 

> 

b" 

7.35  Sec 

b" 

4.14  Sec 

n 

®N  B~  0.63  Sec 


TABLE  2 


Diffusion  Coefficients  for  Binary  Nitrides  (10) 


Nitride 

Dif fusivity 

Temp (K) 

D 

o 

(— ) 
sec 

Q (kJ-Mole 

TIN 

D 

c 

1173-1843 

5.4 

X 

1 

o 

rH 

217.7 

TiN 

D 

c 

1573-1873 

4.4 

X 

<r 

l 

o 

H 

304.7 

TiN 

Dc 

1623-1973 

2.0 

X 

10'3 

376.8 

TiN 

D 

c 

1473-1723 

5.8 

X 

1 

O 

rH 

279.2 

TiN 

DC 

1573-1943 

2.3 

X 

1 

o 

*H 

210.1 

TiN 

DC 

1223-1773 

2.1 

X 

10-7 

199.2 

ZrN 

D 

c 

923-1123 

7.9 

X 

ON 

1 

o 

rH 

150.3 

ZrN 

D 

c 

1623-1973 

6.0 

X 

10-6 

251.2 

ZrN 

D 

c 

1533-1993 

1.7 

X 

io-6 

224.4 

ZrN 

D 

c 

1073-1373 

2.5 

X 

10-6 

185.0 

ZrN 

D 

c 

1473-1773 

3.1 

X 

O 

1 

332.4 

ZrN0.76 

D 

c 

1873-2473 

7.5 

X 

io-5 

327.8 

ZrN0.68 

D 

c 

2823-2973 

3.5 

X 

108 

1025.7 

ZrN0.86 

D 

c 

2673-3023 

7.9 

X 

10 

628.0 

ZrN0.96 

D 

c 

2673-2873 

2.6 

X 

io-2 

406.1 

ZrN0.93 

D 

c 

1273-1473 

4.1 

X 

io-10 

153.2 

ZrN0.75 

D 

c 

1273-1473 

3.0 

X 

io"14 

96.2 

VN 

Dc 

1723 

D = 2 

.5 

x IO-12 

NbN 

D 

1783-2308 

2.1 

468.9 

Table  2 (Con't) 


Nitride 

Dif fusivity 

Temp (K) 

2 

D (— ) 
o sec 

* 

Q (kj-mole 

Cr2N 

D 

c 

1473 

D - 4.2  x 

io'12 

Fe4N 

D 

c 

777 

D - 3.2  x 

io-16 

eFe4N 

D 

c 

853-1003 

4.4  x 

1 

O 

H 

113.4 

ThN 

D 

c 

2173-2673 

2.5  x 

io'2 

416.1 

UN 

°T 

1773-2173 

2.6  x 

00 

• 

O 

H 

230.2 

UN 

dt 

2065 

D ■ 8.1  x 

io*15 

UN 

dt 

1973-2273 

-v-234.4 

UN 

D 

c 

2073-2673 

1.2  x 

io'3 

502.4 

aSi-jN. 

■i  4 

°T 

1473-1683 

1.2  x 

io-16 

233.2 

6-Si3NA 

dt 

1473-1683 

5.8  x 

io2 

777.4 

S13N4 

D 

c 

1763-2023 

io6 

778.7 

As  already  pointed  out,  the  lack  of  a long  lived  nitrogen  radio- 
tracer  has  led  to  difficulties  in  making  nitrogen  transport  measure- 
ments. Chemical  nitrogen  dlffuslvltles  (calculated  by  following  the 
changes  in  a property),  Dc»  have  been  determined  by  exposing  a sample 
to  a fixed  environment  and  following  either  the  kinetics  of  weight 
change,  electrical  conductivity  change,  or  density  change.  By  apply- 
ing the  appropriate  equations  (11-14)  the  diffusivity  can  be  calcu- 
lated assuming  the  mechanism  of  the  process  is  known  exactly.  These 
methods  are  subject  to  considerable  question  due  to  the  influence  of 
grain  boundary  diffusion  contributions  and  the  assumptions  required 
by  the  model. 

Tracer  dif fusivities,  D^,,  can  be  obtained  by  forming  a diffusion 
couple  where  one  of  the  segments  is  enriched  in  nitrogen-15  related 
to  the  other.  Alternatively  a sample  is  placed  in  a gaseous  envi- 
ronment enriched  in  nitrogen-15  and  the  decrement  of  nitrogen-15  in 
the  gas  phase  during  the  exchange  is  determined.  This  method  can 
give  erroneous  results  due  to  exchange  with  the  containment  vessel. 

If  the  sample  is  polycrystalline  a way  of  separating  the  contribution 
of  grain  boundary  diffusion  from  the  volume  diffusion  must  be  de- 
veloped. It  is  also  possible  to  obtain  a tracer  diffusivity  by  mea- 
suring the  nitrogen-15  concentration  versus  depth  profile  after  ex- 
change. The  concentration  versus  depth  profile  can  either  be  ob- 
tained by  mass  spectrometry  or  activation  analysis . To  analyze  a 

sample  by  mass  spectrometry,  the  sample  must  first  be  sectioned  and 

14  15 

converted  to  a gaseous  phase  where  the  N - N ratio  is  measured. 
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This  procedure  has  two  possible  sources  of  error,  which  become  in- 
creasingly important  when  attempting  to  measure  small  penetrations. 
These  are  in  determining  the  exact  amount  of  material  contained  in 
each  section  and  in  determining  the  degree  of  conversion  to  a gas- 
eous phase.  Activation  analysis  can  be  performed  either  by  induced 
radioactivity  or  the  direct  observation  of  a nuclear  reaction.  (The 
direct  observation  of  a nuclear  reaction  technique  is  discussed  in 
Appendix  A.)  In  1967  Condit,  Holt,  and  Hiramel  (15)  suggested  an  acti- 
vation analysis  technique  for  detecting  nitrogen-15  in  solids  using 
15  18 

the  N(a,n)  F reaction.  The  reaction  product  Fluorine-18  is  radio- 
active, which  decays  by  emission  of  a 0.67  MeV  B+  with  a half-life 
of  6732  seconds . The  B+  particles  are  then  detected  by  autoradio- 
graphy. This  technique  was  applied  by  Holt  and  Almassy  (16)  to  study 
nitrogen  diffusion  in  uranium  nitride  in  1969. 

3.2  Nitrogen  Transport  in  Silicon  Nitride 

To  date  there  have  been  few  studies  which  have  allowed  deter- 
minations of  nitrogen  dif fusivities  in  silicon  nitride.  Table  3 (17) 
tabulates  the  activation  energies  reported  to  date  for  nitrogen 
transport  in  silicon  nitride  and  the  method  by  which  it  was  obtained. 
The  large  variations  in  activation  energies  can  possibly  be  attribu- 
ted to  differences  in  the  impurity  levels  in  the  samples. 

Kijima  and  Shirasaki  (17)  obtained  tracer  diffusion  coefficients 
for  both  a-  and  p-Si^N^  by  gas-solid  exchange  techniques.  The 
decrement  in  nitrogen-15  isotope  in  the  gas  phase  was  followed  as  a 
function  of  time  during  the  exchange  and  the  data  reduced  using  the 


TABLE  3 


Activation  Energies  for  Nitrogen  Transport  in  Silicon  Nitride 


Temp  Range 
(K) 

Activation  Energy 
-1 

(kJ  - mole  ) 

Method 

Ref 

1243-1873 

141 .5 

Nitriding  Kinetics 

18 

1473-1683 

233.2 

Isotopic  Exchange 

17 

1533-1683 

418.6-837.3 

Nitriding  Kinetics 

19 

1623-1683 

661.5 

Nitriding  Kinetics 

20 

1473-1683 

777.4 

Isotopic  Exchange 

17 

1763-2023 

778.7 

Thermal  Decomposition 
Kinetics 
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method  of  Carmon  and  Haul.  Their  samples  were  prepared  by  nitriding 

high  purity  silicon  pellets  to  produce  pure  a-Si^N^.  The  6-Si^N^ 

4 

was  produced  by  firing  the  a-Si^N^  at  1653  K for  4.32  x 10  seconds. 
The  two  types  of  Si^N^  are  reported  to  have  similar  densities,  55% 
theoretical,  but  to  have  different  grain  sizes.  The  average  grain 
size  for  the  a-Si^N^  was  50.0nm  and  was  150. Onm  for  the  0-SLjN^. 

The  pellets  were  .rushed  to  form  powders  with  average  particle  sizes 
of  75pm  and  1.0mm  which  were  used  for  the  exchange.  It  was  necessary 
to  substitute  the  average  grain  size  for  the  average  sphere  size  in 
the  equation  for  the  early  stages  of  exchange  between  a sphere  and  a 
well-stirred  solution  of  limited  volume  in  order  to  obtain  agreement 
between  the  dif f usivities  calculated  for  the  two  particle  sizes. 

The  measured  dif f usivities , corrected  for  grain  size,  for  the  temper- 
ature range  1473  to  1673  K are: 


[1.2  x 10  exp  - C 


,2.33  x 


io5j 


RT 


>] 


ID 

sec 


(32) 


for  a-Si^N^  and: 


5 2 

r r rt  ,7.77  X 10  J,,  B 

[5.8  x 10  exp  -( )]  — c 


(33) 


for  6-Si-N. . 
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Kijima  and  Shirasaki  report  that  minimal  experimental  error  was 
introduced  by  reactions  (thermal  decomposition,  oxidation,  reaction 
between  the  sample  and  crucible  and  nitriding  of  the  crucible  or 
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suscepter).  This  was  determined  by  monitoring  the  relative  pressure 

of  nitrogen  against  a standard  of  argon  by  mass  spectrometry  and  also 

checking  for  weight  change  of  the  nitride  powder,  crucible,  and  suc- 

3 

cep tor  during  a 3.6  x 10  second  anneal  at  1623  K.  It  would  seem 
doubtful  that  either  mass  spectrometry  or  weight  change  measurements 
would  be  sensitive  enough  to  detect  the  presence  of  a reaction  during 
an  anneal  of  this  length.  Because  no  reference  is  made  as  to  sample 
size  or  to  the  amount  of  gas  present  the  following  assumptions  will 
be  made:  particle  size  = 75  pm  spheres,  initial  mass  of  Si^N^  « 10  mg, 
and  volume  of  gas  present  * 1 liter.  Assuming  that  l.Onm  of  SiO 


forms  on  each  particle  the  weight  change  would  be  0.1  pg  and  the 
volume  of  nitrogen  would  change  by  2 x 10  "*%.  Since  all  the  particles 
will  not  be  equally  exposed  to  the  nitrogen  these  changes  will  prob- 
ably be  even  smaller. 

It  is  also  emphasized  that  Kijima  and  Shirasaki's  technique 
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allowed  nitrogen  dif fusivities  less  than  10  m - sec  to  be 
determined.  The  diffusion  times  quoted  by  Kijima  and  Shirasaki  are 


of  the  order  of  5 x 10  seconds . This  corresponds  to  a maximum 

penetration  distance  of  2 nm  and  to  0.07  nm  for  their  lowest  diffusi- 
-24  2 -1 

vlty  ('vlO  m - sec  ) which  is  considerably  less  than  a unit  cell 


dimension  of  B-Si^N^.  It  seems  improbable  that  any  technique  would 
allow  the  measurement  of  diffusion  coefficients  corresponding  to 


penetration  distances  on  the  order  of  a unit  cell.  Even  if  it  were 
possible  to  measure  penetration  distances  to  a unit  cell  dimension 
or  even  several  unit  cells,  the  calculated  diffusivity  would  not  be 
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one  representative  of  the  bulk  diffusivity  (22),  but  will  be  charac- 
teristic of  the  sample's  surface.  For  penetrations  this  small  the 
measurement  is  probably  only  representative  of  a phase  boundary  re- 
action. 

Wuensch  and  Vasiles  (23)  proposed  to  study  both  silicon  and 
nitrogen  self-diffusion  in  silicon  nitride  by  ion  microprobe  spectro- 
metry, but  were  unsuccessful  due  to  charge  buildup  on  the  sample  sur- 
face while  sputtering.  They  therefore  microchemically  separated 
the  silicon  and  nitrogen,  and  finally  analyzed  the  separated  com- 
ponents for  isotopic  ratio  by  mass  spectrometric  techniques.  Due  to 
carbon  contamination  during  sectioning,  nitrogen  isotopic  gradients 
were  "not  as  satisfactory  as  hoped"  and  no  self-diffusion  coefficient 
for  nitrogen  was  obtained. 

Brook,  et.al.  (24)  calculated  a process  diffusion  coefficient 

from  densif ication  rates  of  Si^N^  containing  5W/o  MgO  during  hot 
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pressing  at  21  MN-m  . They  assumed  that  the  hot  pressing  model 
proposed  by  Coble  (14)  could  be  used  where  the  rate  of  densif ication 
is  controlled  by  grain  boundary  diffusion.  The  model  has  been  modi- 
fied for  liquid-phase  sintering  to  allow  use  of  the  following  simpli- 
fied rate  equation: 


dp 

dt 


47V6D 


k Ta 


GB  p 
3 *A 


(34) 


where 


dp. 

dt 


is  the  densif ication  rate,  V is  the  volume  transported  by 
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each  atom  of  the  least  mobile  species,  and  PA  is  the  applied  stress. 

The  grain -boundary  thickness,  6,  was  calculated  from  the  grain  size 
and  the  amount  of  MgO  added  (it  was  assumed  all  of  the  MgO  combined 
with  SiO^  to  form  grain  boundaries)  on  the  basis  of  a uniform 
boundary  film.  Using  this  model  they  found  an  activation  energy 
for  nitrogen  grain-boundary  diffusion  of  4.5  x 10^  J-mole  ^ for  the 
temperature  range  1823-2023K  and  6.95  x 10"*  J-mole  * for  the  tempera- 
ture range  1723-1823K.  They  interpret  the  break  at  1823  K as  being 
due  to  the  solidus  of  the  MgSiO^  ~ SiO^  system  (1916K)  with  diffusion 
being  in  the  second  phase  boundary  either  containing  or  not  contain- 
ing a liquid  phase. 

It  would  appear  as  if  the  dif fusivities  calculated  by  this  model 
(10  - 10  1^  m^-sec  are  several  orders  of  magnitude  too  low  con- 

sidering they  represent  diffusion  in  either  a liquid  phase  or  a phase 
that  is  almost  at  its  melting  point.  There  is  no  data  available  for 
nitrogen  diffusion  in  molten  silicates,  however  one  can  assume  it 
should  be  of  approximately  the  same  order  of  magnitude  as  other 
dif fusivities  in  molten  silic  tes.  Oishi,  et.al.  (25)  have  quoted 

dif fusivities  for  oxygen,  calcium,  aluminum,  and  silicon  in  molten 
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silicates  to  be  of  the  order  of  10  m -sec.  The  reason  for  the 
low  diffusivities  can  possibly  be  attributed  to  two  equilibrium 
assumptions  required  to  use  the  model.  The  first  assumption  involves 

% 

! the  grain  size  of  the  hot  pressed  material;  the  diffusivity  is  pro- 

portional to  the  cube  of  the  grain  size.  During  the  hot  pressing 


process  considerable  grain  growth  occurs  as  the  Si^N^  density  in- 
creases. The  model  requires  knowledge  of  the  grain  size  at  the  time 
the  measurement  is  made  and  could  be  rapidly  changing.  The  diffusion 
coefficient  is  inversely  proportional  to  the  grain  boundary  thickness 
which  has  been  calculated  from  the  grain  size  and  the  equilibrium 
phase  diagram  for  the  MgSiO^  - Si02  system.  In  making  this  calcula- 
tion it  is  therefore  necessary  to  know  the  grain  size,  which  is  not 
known  exactly  as  discussed  above,  and  to  assume  all  of  the  MgO  added 
goes  to  the  grain  boundaries  where  it  combines  with  SiC^.  Even  if 
all  of  the  MgO  does  go  to  the  grain  boundaries  the  grain  boundary 
film  is  probably  far  from  being  in  equilibrium. 


SECTION  IV 


Experimental 

4.1  Selection  of  a Technique 

It  was  desired  to  be  able  to  measure  the  self-diffusion  coefficient 
for  nitrogen  in  both  grains  and  grain  boundaries  of  silicon  nitride. 

The  stable  isotope,  nitrogen-15,  was  used  as  a tracer  and  detected  by 
the  direct  observation  of  a nuclear  reaction.  The  nuclear  micro- 
analytical  technique  eliminates  the  need  for  serial  sectioning  which 
would  create  a serious  problem  if  Kijima  and  Shirasaki’s  (17)  measured 
diffusivity  is  correct.  The  technique  also  eliminates  the  problem 
encountered  by  Wuensch  and  Vasiles  (23)  of  efficiently  converting 
each  section  to  a gaseous  phase,  free  of  contaminants,  suitable  for 
measuring  the  nitrogen  isotopic  ratios. 

There  are  three  nuclear  reactions  available  for  the  detection  of 
nitrogen-15,  these  are  given  in  Table  4 along  with  the  Q value  (Q  is 
the  energy  equivalent  to  the  mass  difference  between  the  reactant 
and  product  nuclei)  and  the  possible  interfernces . The  deuteron  in- 
duced reaction  on  nitrogen-15  was  eliminated  due  to  nitrogen-14  inter- 
ference reactions.  Of  the  two  possible  proton  induced  reactions  on 
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nitrogen-15  the  N(p,aQ)  C reaction  was  selected  because  of  the 
high  Q value  and  the  virtually  background  free,  sharp  resonance  at 
1.210  MeV.  The  reaction  also  has  a broad  resonance  which  could  be 
used  for  the  single  spectrum  technique,  see  Figure  3. 


TABLE  4 


Charged  Particle  Induced  Nuclear  Reactions  on  Nitrogen-15 


Reaction 


Q(MeV) 


4.964 


Interferences 


0.531 


7.683 


. * , : 


4.2  Sample  Selection 


Due  to  the  fact  that  large  single  crystals  of  high  purity,  theo- 
retically dense,  dislocation  free  silicon  nitride  are  not  presently 
available,  commercial  grade  silicon  nitride  was  used  in  these  studies. 
It  was  initially  thought  that  CVD  silicon  nitride  would  make  the  best 
samples.  These  samples  would  have  been  thin  films,  therefore  making 
the  selection  of  the  substrate  material  important  to  prevent  reaction 
with  or  contamination  of  the  film.  It  was  determined  that  silicon 
would  make  the  best  substrate,  but  it  was  found  that  deposition  was 
not  uniform  and  there  was  an  adhesion  problem  with  some  of  the  samples. 

Next  attention  was  turned  to  a piece  of  Norton  reaction  bonded 
silicon  nitride.  The  reaction-bonded  sample,  while  having  a high 
purity,  was  fine  grained  and  was  not  uniformly  dense.  Due  to  the  com- 
plex structure  (small  grains,  high  porosity,  and  nonuniform  density), 
the  models  presently  available  and  those  that  could  possibly  be  de- 
rived were  not  able  to  reduce  the  concentration  versus  depth  profile. 
For  this  reacon  the  reaction-bonded  samples  were  not  used. 

TM 

It  was  finally  decided  to  use  a sample  of  fully  dense  NORALIDE 
NC-132  hot  pressed  silicon  nitride  cut  from  billet  number  122072A. 
Electron  microscopy  studies  (26)  showed  the  sample  to  be  composed  of 
a statistical  mixture  of  ^95%  0.5  pm  grains  of  6-silicon  nitride  and 
,v5%  0.5  pm  grains  of  silicon  oxynitride  surrounded  by  a continous 
glassy  phase,  presumably  a silica  rich  magnesium  oxide  glass  (27-29), 
approximately  3 nra  wide  (29).  See  Figure  4 (26).  The  dark  spots  in 
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Microstructure  of  NORALIDE 
nitride,  10,000  x (26). 


Microstructure  of  NORALIDE  NC-132  hot-pressed  silicon 
nitride,  30,000  x (26). 


the  micrographs  are  tungsten  carbide  (26).  A typical  chemical  ana- 
lysis of  the  material  is  given  in  Table  5 (31). 

The  e form  of  silicon  nitride  is  composed  of  SiN^  tetrahedra  shar- 
ing corners,  with  each  nitrogen  corner  being  common  to  three  tetra- 
hedra. The  tetrahedra  are  arranged  in  a hexagonal  structure  as  shown 
in  Figures  5a  and  5b  (32). 

Silicon  oxynitride,  Si^N^O,  is  composed  ofSiN^O  tetrahedra  ar- 
ranged so  that  parallel  sheets  of  silicon  and  nitrogen  atoms  are  joined 
by  Si-O-Si  bonds.  The  crystal  structure  of  silicon  oxynitride  is 
shown  in  Figure  5c  (32). 

4.3  Sample  Preparation  and  Exchange 

The  samples  were  diamond  machined  from  the  billet  to  a nominal 
size  of  7mm  x 10mm  x 2mm.  The  samples  were  then  mechanically  polished 
using  an  alumina  slurry. 

Thermodynamic  calculations  show  that  at  1600  K an  oxygen  partial 
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pressure  less  than  7.0  x 10  atmosphere  must  be  maintained  at  the 
silicon  nitride  surface  in  order  to  prevent  the  thermal  oxidation 
of  the  surface  while  exchanging  with  one  atmosphere  of  nitrogen.  The 
thermodynamic  calculations  considered  are  given  in  Table  6.  In  order 
to  assure  that  this  oxygen  partial  pressure  was  maintained  at  the 
sample  surface,  the  sample  was  placed  in  a tungsten  sample  basket 
and  completely  surrounded  by  powdered  carbon.  A tungsten  sample 
basket  was  used  instead  of  alumina  for  two  reasons.  First  the  tungsten 
will  act  as  an  oxygen  getter  and  second  the  possiblity  of  reactions. 
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There  are  three  reactions  that  were  considered  (reactions  3-5  in 

Table  6) , while  equilibrium  thermodynamics  will  not  allow  enough 

oxygen  to  be  present  when  A^O-j  is  exposed  to  an  atmosphere  of  Nj  at 

1600  K the  errors  in  Gfj^QO  could  allow  the  oxygen  partial  pressure 

to  exceed  7.0  x 10  ^ atmosphere.  The  other  two  reactions  which 

involve  the  chemical  reaction  between  Si„N.  and  A1„0_  were  of  pri- 

3 4 2 3 

mary  importance  when  considering  the  use  of  thin  film  CVD  Si^N^. 

As  can  be  seen  from  reactions  6-9  in  Table  6 as  long  as  P(c02^1600 
< 6.4  x 10  ^ atmosphere  and  p(C0)^qq  < 1.5  x 10  ^ atmosphere  car- 
bon will  serve  as  an  adequate  getter  to  prevent  the  thermal  oxida- 
tion of  the  silicon  nitride  surface.  Mass  spectrometric  analysis  of 
the  ambient  gas  showed  no  mass  44  peak  and  no  increase  in  the  mass  28 

peak  (it  was  possible  to  use  the  mass  28  peak  since  the  ambient  was 
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99.13%  N thus  giving  a low  ^ peak)  verifying  that  carbon  was 
an  adequate  getter. 

Each  sample  was  heated  in  vacuum  to  the  diffusion  temperature  and 
maintained  in  vacuum  until  a pressure  of  10  ^ atmosphere  was  obtained 
Oxygen  free,  research -grade  nitrogen  (0.37%  ^N)  was  then  introduced 
to  the  system  to  obtain  a pressure  of  one  atmosphere.  The  sample  was 
preannealed  in  this  atmosphere  for  a period  of  time  equal  to  or  great 
er  than  twice  the  planned  diffusion  time.  At  the  conclusion  of  the 
preanneal  the  ambient  was  changed  to  a 99.13%  enriched  nitrogen-15 
gas.  At  the  end  of  the  diffusion  time  the  nitrogen-15  was  removed 
via  zeolite  and  the  sample  cooled  to  room  temperature. 


The  temperature  range  investigated  was  from  1506  K to  1723  K 

3 

with  diffusion  times  ranging  from  5.4  x 10  seconds  to  8.64  x 10 
seconds . 


4.4  Activation  Analysis 

Following  the  diffusion  anneal  the  samples  were  mounted  in 

aluminum  sample  holders,  which  were  mounted  on  a movable  target 

holder  in  the  scattering  chamber  capable  of  holding  five  samples  at 

a time.  The  samples  were  positioned  so  that  the  incident  proton 

beam  made  an  angle  of  7.5°  with  the  sample  normal.  (This  geometry 

was  chosen  to  prevent  channelling  when  analyzing  single  crystals). 

The  energetic  alpha  particles  were  observed  at  a laboratory  scatter- 
o 

ing  angle  of  165  ; after  they  passed  through  a magnetic  spectrometer, 
a momentum  analyzer  (designed  by  Lindstrom  and  Hever  (34)),  by  a 
conventional  silicon  surface  barrier  detector.  The  magnetic  spec- 
trometer eliminates  the  need  for  mylar  to  separate  the  alphas  from 
the  backscattered  protons,  which  is  an  advantage  since  one  does  not 
have  to  worry  about  the  uniformity  or  stopping  power  of  the  mylar 
film.  A schematic  of  the  scattering  chamber  and  magnetic  spectrometer 
is  shown  in  Figure  6. 

Each  sample  was  bombarded  with  a monoenergetic  beam  of  protons 
and  the  resulting  alpha  spectrum  recorded.  By  systematically  vary- 
ing the  incident  proton  energy  upward  from  1.205  MeV  the  1.210  MeV 
resonance  could  be  moved  deeper  into  the  sample.  Figure  7 shows  the 
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mean  depth  associated  with  each  incident  proton  energy.  The  mean 
depth  associated  with  each  incident  proton  energy  was  calculated 
using  the  stopping  powersgiven  by  Williamson,  et.al.  (35).  These 
values  are  accurate  to  within  better  than  10%  (oxide  ceramics  have 
shown  agreement  to  within  5%)  which  produces  a negligible  error  in 
the  final  results.  The  depth  probed  is  determined  by  the  width  of 
the  resonance.  For  incident  proton  energies  less  than  1.223  MeV 
the  depth  probed  is  equal  to  twice  the  mean  depth,  because  the 
resonance  is  not  totally  within  the  sample.  For  energies  greater 
than  1.223  MeV  the  depth  probed  is  0.4  ym  which  is  the  depth  equiva- 
lent to  the  resonant  width. 

Figure  8 shews  the  experimental  alpha  specta  associated  with  the 
1.210  MeV  resonance  in  an  undiffused  (0.37%  ^N)  Si^N^  sample  at 
various  incident  proton  energies.  Although  the  shape  of  the  spectra 
at  the  various  energies  is  different,  due  to  alpha 
straggling,  the  total  area  of  each  peak  is  statistically  the  same, 
4360  ± 67  alphas  per  10,000  yC. 


incident  proton  energies. 


SECTION  V 


RESULTS  AND  DISCUSSION 

5.1  Data  Reduction 

The  number  of  alpha  particles  observed  was  plotted  as  a function 

of  alpha  energy  for  each  incident  energy  and  the  area  (total  number 

of  counts)  associated  with  the  1.210  MeV  resonance  determined.  This 

was  accompl  shed  by  estimating  a background  (the  background  is  due 

to  other  resonances  and  interferences  between  resonances,  see  Section 

B.4)  and  subtracting  it  from  the  number  of  counts  in  the  resonance 

region.  The  background  was  calculated  by  fitting  the  average  number 

of  observed  alphas  over  the  ten  channels  above  the  resonance  region 

and  the  ten  channels  below  the  resonance  regon  to  an  appropriate  curve 

(which  in  most  cases  was  a straight  line).  Figures  9 and  10  are 

representative  of  the  alpha  spectra  obtained  and  the  background 

estimated.  Figure  9 relates  to  a sample  exchanged  with  nitrogen-15 
14 

for  8.64  x 10  seconds  at  1507.5  K.  The  incident  proton  energy  was 

3 

1.300  MeV  and  the  total  integrated  beam  current  was  2 x 10  yC.  The 

area  associated  with  the  1.210  MeV  resonance  was  1357  t 37  counts. 

Figure  10  shows  the  alpha  spectrum  obtained  from  an  undlffused  sample 

(0.372  ^N)  which  was  used  for  normalization  of  the  data.  The  proton 

activation  analysis  parameters  were  the  same  as  in  Figure  9.  The  area 

associated  with  the  1.210  MeV  resonance  was  664  ± 26  counts. 

At  each  energy  the  ratio  of  the  areas  of  the  resonance  region  of 

the  diffused  sample  to  the  undiffused  sample  was  calculated  and  multi- 
-3 

plied  by  3.7  x 10  , the  natural  abundance  of  nitrogen-15,  to 


Observed  alpha  spectrum  for  a sample  exchanged  with  nitrogen 
seconds  at  1507.5  K. 


Figure  10:  Observed  alpha  spectrum  for  an  undiffused 


determine  the  nominal  average  fractional  concentrattoo._Q.f  Jj,ttrQfl3P-l5 , in 
each  depth  Interval.  Figures  11-14  show  the  nominal  average  frac- 
tional concentration  versus  "mean  depth"  profiles  for  representative 
samples.  The  "mean  depth"  associated  with  each  Incident  proton  energy 
was  obtained  from  Figure  7. 

In  order  to  obtain  the  true  average  fractional  concentration, 

C + C , from  the  nominal  average  fractional  concentration  it  is 
o 

necessary  to  look  at  the  variation  in  the  concentration  profile  over 
the  depth  probed  by  the  resonance.  Since  the  yield  from  the  resonance 
region  is  proportional  to  the  Integral  of  the  product  of  the  concen- 
tration and  the  differential  cross  section  over  the  depth  interval 
associated  with  the  depth  of  the  resonance  (see  Equation  38),  the 
nominal  average  fractional  concentration  and  true  average  fractional 
concentration  will  be  exactly  equal  if  the  concentration  profile  in 
the  interval  is  constant.  If  the  concentration  profile  varies  slowly 
in  a uniform  manner  over  the  depth  interval  then  the  nominal  average 
fractional  concentration  will  be  a reasonable  estimate  of  the  true 
average  fractional  concentration.  This  is  the  case  for  the  points 
corresponding  to  "mean  depths"  greater  than  0.25  pm  in  Figures  11-14. 
However,  if  the  concentration  profile  varies  rapidly  or  in  a nonuniform 
manner  across  the  depth  interval  the  nominal  average  fractional  concen- 
tration and  the  true  average  fractional  concentration  can  be  drastic- 
ally different.  This  is  the  case  for  the  first  four  points  in 
Figures  11-14.  The  analysis  of  the  concentration  profile  in  this 
region  will  be  treated  later. 
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Depth  (pm) 

Figure  12:  Nominal  average  fractional  concentration  veraus  depth  profile  for  Sample 


Depth  (/*m) 

Figure  13:  Nominal  average  fractional  concentration  versus  depth  profile  for  Sample 


According  to  Levine  and  MacCallum  (7),  a plot  of  In  C versus  x , 
see  Equation  30,  should  be  linear  at  deep  penetrations.  After  elimi- 
nating the  first  four  points  due  to  the  rapidly  varying  concentration 
in  the  depth  intervals  the  remaining  data  was  plotted  as  In  C versus 
x6/5  and  fit  to  a straight  line  by  a linear  least  squares  procedure 
(36).  To  determine  how  close  to  the  surface  to  extend  the  fit  the 
best  value  of  the  linear  correlation  coefficient  (36)  was  used. 
Figures  15-18  show  representative  fits.  As  can  be  seen  typically  all 
of  the  point  could  be  fit  by  the  line. 

In  order  to  check  the  validity  of  the  model  the  product 


(^6/5)  was  calculated  for  each  sample.  According  to 

3x 

Equation  30  this  product  should  be  a constant  at  any  temperature. 

Table  7 gives  the  experimental  slopes,  times,  and  products  for  the 

two  temperatures,  1565.5  K and  1683.0  K,  where  more  than  one  sample 

SlnC 


was  exchanged.  (Note  that  the  uncertainity  in  the  slope. 


8x 


6/5 


is  typically  less  than  12%. ) As  can  be  seen  from  Table  7 fair  agree- 
ment was  obtained  for  all  samples  with  the  maximum  variation  in  the 
product  being  a factor  of  6.7.  The  lack  of  any  systematic  trend  in 
the  products  with  time  tends  to  indicate  something  inherent  in  the 
samples  themselves  (e.g.  microstructure)  rather  than  an  inadequacy  in 
the  model. 

If  one  takes  the  logarithm  of  each  side  of  Equation  30  it  can  be 

shown  that  ln[  (~^vy ) 2]  should  be  proportional  to  T \ 

3x  ° 


TABLE  7 

Experimental  Parameters  for  Grain  Boundary  Fits 


Sample 

Temperature 

(K) 

,31nC  . 
3x6/5) 

Time 

(Sec) 

<f&5/V 

2 

1568.0 

-2.072 

± 0.212 

4 

0.54  x 10 

(4.04  ± 

0.69)x 

io-3 

3 

1563.0 

-1.387 

± 0.090 

2.16  x 104 

(3.94  ± 

0.43)x 

10~3 

4 

1565.5 

-0.292 

± 0.032 

8.64  x 10  4 

(2.65  ± 

0.48)x 

io-2 

5 

1681.0 

-0.123 

± 0.009 

4 

0.54  x 10 

(4.47  ± 

0.55)x 

io-1 

6 

1685.0 

-0.188 

± 0.021 

2.16  x 104 

(1.10  ± 

0.21)x 

10_1 

This  being  the  case,  a graph  of  this  function  versus  T should  be  a 
straight  line,  however  as  can  be  seen  in  Figure  19, there  is  a fair 
amount  of  scatter  in  the  points.  This  scatter  can  be  attributed  to 
either  nonuniformity  in  the  samples,  or  the  inability  of  the  model 
to  adequately  treat  the  case  of  intersecting  grain  boundaries  or 
possibly  the  presence  of  a phase  boundary  reaction. 

Because  of  the  rapidly  varying  concentration  of  nitrogen-15  at 


the  sample  surface  the  number  of  counts  obtained  in  the  resonance 
region  for  the  first  four  depth  intervals  may  be  almost  entirely  due 
to  this  surface  concentration.  For  this  reason  the  volume  diffusivity 
was  estimated  using  the  total  number  of  counts  obtained  in  the 
resonance  region.  The  use  of  only  the  total  number  of  counts  to 
estimate  Dv  requires  three  assumptions  to  be  made.  First  it  is 
necessary  to  decide  if  a phase  boundary  reaction  plays  any  part  in  the 
profiles  obtained. Because  of  the  inability  to  determine  uniquely  both  a 
value  of  D and  of  K from  the  data  available.it  was  decided  to  assume 

VS 

there  was  no  phase  boundary  reaction  thus  giving  the  lowest  possible 
value  of  Dy.  This  assumption  was  also  made  by  Kijima  and  Shiraslki 
in  their  study. 

It  was  also  necessary  to  make  some  assumption  regarding  the  contri- 
bution to  the  total  number  of  counts  from  nitrogen-15  in  the  grain 
boundaries.  It  has  been  assumed  for  the  purpose  of  calculation  that 
the  contribution  is  negligible, which  is  justified  by  the  fact  that 


I d 


the  apparent  surface  contribution  due  to  nitrogen-15  in  the  grain 
boundaries  (see  Figures  15-18)  is  approximately  1%  while  it  has  been 
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assumed  that  the  surface  concentration  is  100%  (no  phase  boundary 
reaction).  It  should  also  be  noted  that  the  grain  boundaries  account 
for  less  than  2%  of  the  total  sample  volume. 

The  final  assumption  made  is  that  the  entire  complementary  error 
function  is  contained  in  the  depth  interval  being  considered.  This 
assumption  can  only  be  justified  by  the  final  result. 

The  calculation  was  made  by  using  the  relationship  for  the  yield 
of  a nuclear  reaction  given  by: 

N (x  ) ■ N C (x . ) jo  (x.  )Ax  .A£2  (38) 

o 1 1 XX 


Here  N (x.)  is  the  number  of  counts  associated  with  mean  depth  x . ; N 
o l l 

is  the  number  density  of  nitrogen  atoms  in  the  sample;  C(x^)  is  the 
fractional  concentration  of  nitrogen-15  associated  with  mean  depth 
x±;  j is  the  integrated  flux  of  protons;  o(x^)  is  the  differential 
cross  section  associated  with  mean  depth  x^;  Ax^  is  the  depth  inter- 
val associated  with  mean  depth  x_^ ; and  Af2  is  the  solid  angle  of 
observation.  To  obtain  the  total  number  of  counts  in  the  depth  in- 
terval (x^  »x  ) Equation  38  can  be  rewritten  in  integral  form  as: 


Nq  - Nj  Afi 


C(x)a (x)dx 


(39) 


x„ 

i 
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If  the  differential  cross  section  does  not  drastically  change  over 


the  depth  interval  (x^.x^)  it  can  be  replaced  by  an  average  cross 
section  value,  o(x)  and  Equation  39  rewritten  as: 


N m NjAOo(x) 
o 


C (x)dx 


(40) 


4 


in 

fc/l— — 
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For  Si^N^  N is  equal  to  5.5  x 10  nitrogens-m  . The  total  charge 
collected  was  2 mC  which  corresponds  to  an  integrated  proton  flux 
of  1.25  x 1016  protons  and  the  solid  angle  of  observation  was 
7.2  x 10  ^ sr.  Initially  it  was  assumed  that  the  entire  complementary 
error  function  was  contained  within  the  interval  (x^.x^)  - (0,20nm) 
which  is  the  interval  over  which  a proton  loses  1 keV  of  energy  In 
Si^N^.  The  integral  in  Equation  40  can  then  be  replaced  by  (37): 


20nm 


C(x)dx  = 

erfc  j,  dx  = 

J J 

2 (D  tr 

erf  c 


dx 


2(Dvt) 


1.1284(Dvt) 


(41) 


This  substitution  can  be  made  if  the  surface  concentration  is  unity 
and  the  entire  complementary  error  function  is  contained  in  the  in- 
terval (0,20nm).  Equation  40  can  now  be  rewritten  as: 


Nq  - (5.58  x 1040)  o(x)(Dvt)is 
m 


(42) 
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The  values  of  cr(x)  used  (see  Figure  3)  were  114  mb,  157  mb,  and 
127  mb  for  incident  proton  energies  of  1.205  MeV,  1.210  MeV,  and 
1.215  MeV  respectively.  The  appropriate  number  of  counts  (observed- 
background-undif fused)  and  time  were  substituted  in  Equation  42  and 
the  Equation  solved  for  Dv-  If  it  was  found  that  the  value  of 
2(Dvt)  exceeded  20nm,  as  was  the  case  for  samples  4 and  6,  it  was 
necessary  to  increase  the  depth  interval  and  hence  use  a slightly 
different  differential  cross  section.  This  calculation  was  made  for 
the  first  three  points  of  each  sample  and  the  results  tabulated  in 
Table  8.  As  can  be  seen  reasonable  agreement  was  obtained  between 
the  values  of  Dv  calcualted  from  the  resonance  yield  in  each  of  the 
first  three  depth  intervals  considering  the  fact  that  the  cross 
section  is  known  to  ±5%  at  the  resonant  energy.  The  fact  that  the 
calculated  value  of  Dv  increases  with  increasing  time  at  a given 
temperature  (Samples  2,  3,  & 4 and  5 & 6)  is  indicative  of  either 
the  presence  of  a phase  boundary  reaction  or  the  failure  to  account 
for  the  nitrogen-15  present  in  the  grain  boundaries.  It  can  also  be 
seen  that  the  value  of  calculated  when  the  incident  proton  energy 
was  1.210  MeV  is  usually  the  lowest  value  of  Dv-  This  can  be  at- 
tributed to  the  fact  that  the  differential  cross  section  is  a maxi- 
mum at  this  point  and  any  deviation  in  machine  energy  will  cause 
the  differential  cross  section  to  decrease,  thus  decreasing  the 
alpha  yield. 

Figure  20  shows  a comparison  of  the  number  of  counts  observed  in 
each  of  the  first  three  depth  intervals  and  the  number  of  counts 
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Figure  20:  Comparison  of  the  observed  number  of  counts  and  predicted  number  of  counts  based 
on  the  average  value  of  D given  in  Table  8 for  (a)  Sample  7 and  (b)  Sample  5. 
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observed  in  each  of  the  first  three  depth  intervals  and  the  number  of 
counts  that  would  be  predicted  in  the  interval  assuming  the  average 
diffusivity  given  in  Table  8.  Figure  20a  is  for  the  sample  with  the 
smallest  spread  in  dif fusivities  (Sample  7)  and  Figure  20b  for  that 
with  the  largest  spread  in  calculated  dif fusivities  (Sample  5) . 


The  temperature  dependance  of  Dy  is  shown  in  Figure  21  and  given 


by: 


2 ..  +4.90  x 10-9 

Dv(iic)  = 2‘02  x 10  (-2.0!  x 10~ll)exp 


(3.17  t 0 . 74)x  10' 
RT 

(43) 


The  relationship  found  by  Levine  and  MacCallum  (7),  Equation  30, 
was  used  to  obtain  the  value  of  D 6 for  each  sample.  The  value  of 

CjD 

D substituted  in  Equation  30  was  the  average  value  given  in  Table  8. 
V 9lnC 


Table  9 is  a tabulation  of  the  value  of  . , 6/5  used,  the  values 

a tyt  ) 

of  the  dimensionless  parameters  y , t,  and  yx  , and  the  calculated 
value  of  Dgb  (with  6 = 3.0  nm) . The  value  of  the  grain  boundary 
thickness  used,  6 = 3.0  nm,  was  measured  by  electron  microscopy  (26) 
and  is  in  agreement  with  the  value  given  in  the  literature, 

3.0  ± 1.0  nm  (28). 

The  temperature  dependence  of  D^g  is  shown  in  Figure  22  and  given 
by: 

_ ,m2  , . „„  ,+2.34  x 106x  (5.09  ± 1.77)x  105 

DGbW  ' 4’09  ^ )**P  ‘ 


-4.09 


RT 


(44) 


Figure  21:  Temperature  dependence  of  D 
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5.2  Comparison  with  the  work  of  Kilim  a and  Shirasaki  (17) 


Figure  23  shows  a graphical  comparison  of  Kijima  and  Shirasaki's 

results  for  nitrogen  diffusion  in  6~Si and  the  results  obtained  in 

this  study.  As  can  be  seen  the  value  of  from  this  study 

is  on  the  average  1^  orders  of  magnitude  greater  than  that  calculated 

by  Kijima  and  Shirasaki  while  the  activation  energy  is  approximately 
2 

■j  that  calculated  by  Kijima  and  Shirasaki. 

The  discrepancy  in  the  magnitude  of  the  measured  dif fusivities 
can  be  attributed  to  the  difference  in  starting  materials.  Kijima 
and  Shirasaki  started  with  high  purity,  55Z  theoretically  dense,  reac- 
tion sintered  g-Si^N^  powders.  The  particle  sizes  in  the  powders 
were  75.0  pm  and  1.0  mm  with  average  grain  sizes  of  150  nm.  The 

present  study  used  a slab  of  commercial  grade  hot  pressed  Si^N^, 

TM 

NORALIDE  NC-132  , which  contained  a statistical  mixture  of  95% 

(5  -Si^N^  grains  and  5%  si2N2^  Srains  with  an  average  grain  size  of 
0.5  pm.  The  grains  were  surrounded  by  a 3.0  nm  glassy  phase. 

Gauckler  and  Lucas  (38)  have  found  that  low  concentrations  of  metal 
oxides  will  form  single-phase  solid  solutions  with  g-Si^N^  without 
the  production  of  extrinsic  defects  by  maintaining  a metal  to  non- 
metal  ratio  of  3 to  4.  They  argue  that  the  high  degree  of  covalency 
of  the  Si-N  bond  is  responsible  for  the  solid  solution  formation  and 
the  lattice  not  forming  extrinsic  defects  upon  the  addition  of  metal 
oxides.  (The  justification  for  this  is  not  obvious.)  It  has  been 
found  that  the  g-Si^N^  lattice  expands  when  metal  atoms  larger  than 
silicon  dissolve  in  the  lattice  and  contracts  when  metal  atoms 
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Figure  23:  Comparison  of  D 's  obtained  in  thi9  study  and  by  Kijima 


atoms  smaller  than  silicon  dissolve  In  the  lattice.  It  would  also  be 
possible  that  the  addition  of  MgO  will  create  either  magnesium  inter- 
stitials or  nitrogen  vacancies.  It  is  therefore  reasonable  to  argue 
that  based  on  either  Gauckler  and  Lucas’  argument  (Mg  covalent  radii  = 
0.14  nm  and  Si  covalent  radii  ■ 0.12  nm)  or  the  formation  of  nitrogen 
vacancies  that  the  rate  of  nitrogen  diffusion  will  be  enhanced  by  the 
addition  of  MgO. 

5.3  Comparison  with  the  Work  of  Brook,  et.al.  (24). 

A comparison  of  the  grain  boundary  diffusion  coefficients  calcula- 
ted by  Brook,  et.al.  and  those  measured  in  this  study  is  shown  in  Figure 
24.  (Also  shown  is  the  value  obtained  if  Kijima  and  Shirasaki's 
value  of  Dv  is  used  in  the  calculation).  The  grain  boundary  diffusivity 
measured  in  this  study  is  approximately  1 ^ orders  of  magnitude  larger 
than  that  measured  by  Brook,  et.al.  while  the  activation  energy  is  in 


< 
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experimental  agreement  with  that  of  Brook,  et.al.  Part  of  this  dis- 
crepancy can  be  attributed  to  the  materials  used.  Both  studies  used 
hot  pressed  silicon  nitride,  however  the  amount  of  MgO  used  as  a hot 
pressing  additive  was  different.  Brook,  et  al's  samples  contained 
5 w/0  MgO  and  those  in  this  study  contained  0.5  w/0  MgO. 

Most  of  the  discrepancy  can  be  attributed  to  the  models  used  to 
analyze  the  raw  data.  Brook,  et.al.  used  a modified  version  of  Coble's 
model  (13)  for  densif ication  during  hot  pressing  which  requires 
several  assumptions  (see  Section  3.2)  and  to  the  author’s  knowledge 
results  using  the  model  have  never  been  compared  with  other  grain 
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boundary  diffusion  data.  Of  the  assumptions  made  several  are  subject 
to  considerable  question.  First,  is  nitrogen  grain  boundary  diffusion 
the  rate  controlling  process  in  densif ication  and  if  so,  is  the  modi- 
fied Coble  model  valid?  The  model  used  by  Brook,  et.  al.  also  requires 
a knowledge  of  the  grain  size  and  the  grain  boundary  thickness  at  the 
time  of  the  measurement.  Brook,  et.  al.  do  not  state  how  they  esti- 
mated their  grain  size  especially  since  grain  growth  is  occurring 
during  the  process.  It  was  then  assumed  that  by  a knowledge  of  grain 
size,  the  amount  of  MgO  added,  and  the  MgSiO^-SiC^  equilibrium  phase 
diagram  the  grain  boundary  thickness  could  be  calculated.  As  can  be 
seen  from  Equation  34  the  grain  boundary  diffusivity  is  proportional 
to  the  cube  of  the  grain  size  divided  by  the  grain  boundary  thick- 
ness. Assuming  a constant  grain  boundary  volume,  it  can  be  shown 
that  if  the  grain  size  is  underestimated  by  a factor  of  4 the  grain 
boundary  thickness  will  also  be  underestimated  by  a factor  of  4 which 
will  underestimate  the  grain  broundary  diffusion  coefficient  by  a 
factor  of  16. 

In  the  present  sutdy  it  is  not  known  what  effect  a phase  boundary 
reaction  at  either  the  grain  boundary-surface  or  grain  boundary-grain 
interface  will  have  on  the  value  of  D„„.  It  is  also  not  known  how 
the  presence  of  5Z  silicon  oxynitride  grains  along  with  the  8-Si^N^ 
grains  will  have  on  D__. 

Ufi 

5.4  Conclusions 

This  study  has  shown  that  the  direct  observation  of  the 

15  12 

N(p,a  ) C reaction  can  be  used  to  measure  nitrogen-15  concentration 
0 
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versus  depth  profiles  with  submicron  resolution  and  without  the  need 
for  serial  sectioning.  The  technique  has  been  applied  to  the  study 
of  nitrogen  diffusion  in  NORALIDE^  NC-132  hot  pressed  silicon  nitride 
and  values  of  both  the  volume  and  grain  boundary  dif fusivities  esti- 
mated. 
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SECTION  VI 


FUTURE  WORK 

The  obvious  extention  of  this  work  is  to  investigate  nitrogen 
diffusion  in  both  pure  a-  and  pure  6-Si^N^  over  as  large  a tempera- 
ture range  as  possible.  Once  a diffusion  coefficient  is  obtained  for 
the  pure  forms  of  silicon  nitride  the  effect  of  impurities  should  be 
investigated,  especially  those  added  during  hot  pressing.  If  possible 
a glass  of  the  same  composition  as  that  found  in  the  grain  boundaries 
should  be  produced  to  allow  the  study  of  nitrogen  diffusion  in  a 
material  with  the  same  composition  as  the  grain  boundary. 

In  order  to  eliminate  the  question  of  a phase  boundary  reaction 
it  would  be  advantageous  to  develop  a method  of  producing  layered 

samples  of  theoretically  dense  silicon  nitride.  By  the  samples  being 

14  15 

composed  of  a layer  of  Si^  followed  by  a layer  of  Si^  the 

possibility  of  a phase  boundary  reaction  would  be  eliminated. 

The  proton  activation  technique  should  also  be  applied  to  the 
problem  of  the  nitriding  of  silicon.  By  a method  similar  to  the  one 
used  in  this  laboratory  to  study  the  mechanism  of  the  thermal  oxida- 
tion of  silicon  (37)  it  would  be  possible  to  determine  the  mechanism 
of  the  nitridation  of  silicon. 

A deconvolution  program  would  be  valuable  in  allowing  utilization 
of  the  1.028  MeV  resonance  to  obtain  concentration  versus  depth  pro- 
files by  the  single  spectrum  technique. 

The  scattering  chamber  should  also  be  modified  to  allow  the  inci- 

dent beam  to  strike  the  sample  surface  at  a shallow  angle  (38).  This 


would  be  especially  useful  for  shallow  penetrations , since  the  depth 
resolution  could  be  reduced  to  2-3  nm. 


Finally,  the  proton  activation  technique  should  be  applied  to 
other  nitrogen  containing  systems  and  also  to  some  non-nitrogen 
systems  which  have  a high  nitrogen  solubility  at  their  forming  tem- 
peratures . 
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emission  for  qualitative  and  quantitative  spectrochemical  analysis 
was  in  1913  when  Moseley  (42)  showed  there  was  a relationship  be- 
tween atomic  number  and  the  wavelength  of  the  emitted  X-ray.  How- 
ever, for  a long  time  analytical  x-ray  emission  analysis  was  mainly 
by  x-ray  fluorescence  even  though  it  was  known  that  there  was  a high 
cross  section  for  x-ray  production  due  to  charged  particle  bombard- 
m'nt.  In  1970  Johansson,  et.al.  (43)  showed  experimentally  that 
high  sensitivity,  multielement  analysis  was  possible  by  proton  in- 
duced x-ray  emission. 

Theory  Valkovic  (44)  has  probably  provided  the  best  treatment 
of  the  fundamentals  of  P1XE,  the  basis  of  which  is  given  here.  If 
by  some  means  a vacancy  is  created  in  one  of  the  inner  electron 
shells  it  will  be  immediately  filled  ( 10  ^sec)  by  an  electron  from 
one  of  the  outer  shells  with  the  simultaneous  emission  of  a charac- 
teristic x-ray.  The  x-rays  emitted  when  an  electron  fills  a vacancy 
in  the  K-shell,  n=l,  result  in  the  x-ray  lines  of  the  K-series.  If 
the  transition  is  from  the  L-shell,  n=2,  the  associated  x-ray  is 

known  as  K and  if  the  transition  is  from  the  M-shell,  n=3,  the 
a 

associated  x-ray  is  known  as  K . The  x-ray  energy  increases,  wave- 
rs 

length  decreases,  from  the  M-series  to  the  L-series  to  the  K series 
due  to  the  energy  difference  between  levels.  As  the  atomic  number 
increases,  the  energy  of  the  same  x-ray  spectral  line  also  increases 


due  to  the  increase  in  nuclear  charge,  the  electron  binding  energy 
and  the  increase  in  energy  difference  between  the  orbits.  This 
relationship  is  shown  in  Figure  A.l.  It  is  also  of  importance  to 
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note  this  because  as  the  x-ray  spectra  originates  in  the  inner  orbi- 
tals it  is  almost  completely  independent  of  the  chemical  state  of 
the  atom. 


A continuous,  bremsstrahlung , spectrum  of  x-rays  will  exist  in 
addition  to  the  line  spectrum  described  above.  The  continuous  spec- 
trum is  due  to  electrons  and  other  charged  particles  slowing  down 
in  a step  by  step  fashion.  The  highest  energy  of  the  continuum  will 
correspond  to  deceleration  of  the  highest  energy  particle  to  zero 
velocity. 

Sensitivity  It  is  of  importance  to  examine  how  various  experi- 
mental parameters  effect  the  sensitivity,  minimum  detectable  con- 
centration, keeping  in  mind  that  the  sensitivity  refers  to  the  sam- 
ple being  bombarded.  Sample  preparation  techniques  may  allow  the 
sensitivity  to  be  enhanced.  Once  one  knows  the  minimum  detectable 
concentration,  it  is  possible  to  calculate  the  minimum  detectable 
absolute  amount  for  various  elements  knowing  the  mass  of  the  material 
Irradiated.  Assuming  the  incident  particles  are  protons  in  the  MeV 


{ 

.4 


range,  a sufficient  number  of  x-ray  counts  can  be  obtained  for  masses 
-19 

as  small  as  10  kg.  However,  the  trace  element  to  be  measured  is 
always  contained  in  a matrix,  and  the  counts  associated  with  the 
trace  element  must  be  separated  from  these  due  to  the  bremsstrahlung 
background  and  the  matrix. 

Johansson  and  Johansson  (45)  have  shown  that  the  sensitivity  of 
the  experiment  varies  as  the  square  root  of  the  detector  resolution 
and  inversely  as  the  square  root  of  the  product  of  the  solid  angle 
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subtended  by  the  detector,  the  integrated  flux  of  incident  particles 
and  the  target  thickness.  In  any  experiment  one  is  free  to  choose 
the  bombarding  energy,  the  particle  type,  the  beam  current,  the 
analysis  time,  etc.  in  order  to  optimize  the  sensitivity  for  the 
elements  of  Interest.  For  practical  reasons  the  analysis  time  is 
usually  minimized  and  the  beam  current  is  limited  by  the  nature  of 
the  target,  therefore  one  usually  selects  an  incident  ?nergy  and 
particle  type  to  give  the  best  overall  sensitivity.  Figure  A. 2 (45) 
shows  the  minimum  detectable  fractional  concentration  for  incident 
protons  of  incident  energy  1 to  10  MeV  as  a function  of  atomic  num- 
ber. The  parameters  used  in  Figure  A. 2 relate  to  a carbon  matrix  of 

-3  -2  -2 

thickness  10  mg-mm  , a detector  subtending  an  angle  at  3.8  x 10 

sr  with  a resolution  of  165  eV  and  a total  integrated  proton  flux 

of  10  yC. 

Experimental . A typical  scattering  chamber  for  PIXE  analysis  is 

shown  in  Figure  A. 3 (46).  A monoenergetic  particle  beam  is  passed 

through  a diffuser  (usually  havar  foil)  and  then  through  a colli- 

minator  in  order  to  produce  a homogenous  beam.  The  beam  then  strikes 

o 

the  sample  at  an  angle  of  45  to  the  sample  normal  and  the  character- 
istic x-rays  are  observed  at  a laboratory  angle  of  90°. 

Sample  preparation  produces  the  most  difficulties  in  PIXE  ana- 
lysis. Thin  targets  are  desired  in  order  to  gain  sensitivity,  how- 
ever certain  materials  such  as  biological  tissue  and  metallurgical 
samples  are  difficult  to  obtain  in  thin  sections.  The  problem  of 
using  thick  targets  has  been  treated  by  Patnaik  and  Dhere  (47).  In 
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most  cases  the  material  to  be  analyzed  Is  deposited  on  a suitable 
backing,  such  as  carbon,  collodion,  formuar,  kapton,  mylar,  poly- 
styrene, millipore,  nuclopore,  etc. 

It  is  sometimes  desirable  to  analyze  liquid  samples  or  samples 
which  deteriorate  in  vacuum.  In  these  cases  external  beams  can  be 
used.  The  particle  beam  is  passed  through  a thin  window,  usually 
mylar  (48)  nickel  (48),  or  beryllium  (49),  which  sacrifices  sensi- 
tivity due  to  increased  bremsstrahlung,  but  simplifies  sample  pre- 
paration. 

Applications.  The  two  primary  uses  of  PIXE  have  been  in  the  areas 

of  environmental  studies  (50-56)  and  biological  systems  (56-63) 

which  include  plant  material  (56-60),  animal  tissue  and  fluids 

(57-63),  and  medical  uses  (6-63).  The  application  to  materials 

systems  is  presently  in  the  developmental  stages  due  to  the  problems 

associated  with  thick  samples.  Research  is  presently  being  carried 

out  which  will  allow  concentration  versus  depth  measurements  to  be 

made  (64-66)  by  one  of  three  methods  1)  variation  of  incident  ion 

energy;  2)  variation  of  the  angle  between  the  incident  beam  and  the 

target,  and  3)  use  of  the  ratio  between  the  K„  and  K x-ray  intensi- 

p a 

ties. 

A. 3 Rutherford  Backscattering 

Introduction.  Rutherford  Backscattering,  RBS,  is  the  oldest  of 
the  nuclear  microanalytical  techniques.  The  first  application  of 
RBS  to  materials  analysis  was  by  Geiger  and  Marsden  (67)  in  1909 
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and  was  subsequently  explained  In  a paper  by  Rutherford  (68)  in 
1911.  The  technique  was  not  well  publicized  until  the  1960’s 
when  a 5 MeV  alpha  source  was  part  of  the  Surveyor  5 payload  (69,70). 

The  backscattered  alphas  were  used  to  provide  Information  about  the 
elemental  composition  of  the  moon. 

RBS  is  a simple  technique  which  is  absolute,  it  requires  no  addi- 
tional calibration,  and  will  give  reliable,  quantitative  elemental 
depth  distributions  provided  the  sample  surface  is  smooth  and  the 
elelments  to  be  considered  do  not  have  similar  masses. 

Theory.  Chu  (71)  has  divided  RBS  into  three  basic  concepts:  1) 
the  kinematic  factor  (qualitative  analysis),  2)  the  differential 
scattering  cross  section  (quantitative  analysis),  and  3)  the  energy 
loss  (depth  analysis).  Each  of  these  concepts  will  be  discussed 
below. 

0 

The  kinematic  factor,  K,  is  defined  as  the  ratio  of  the  scattered 
particle's  energy,  Eg,  to  its  energy  before  scattering,  E.  The 
kinematic  factor  can  be  simply  calculated  from  the  laws  of  conser- 
vation of  mass-energy  and  conservation  of  linear  momentum.  If  the 
mass  of  the  incident  particle,  ®,  the  mass  of  the  target,  and 
the  laboratory  scattering  angle,  0,  are  known  and  the  process  is  as- 


The  differencial  scattering  cross  section. 


do  (6) 


, arises  from 


a pure  coulomb  interaction  between  the  incident  particle  and  the 
target  and  is  a classical  physics  problem  (72).  The  differential 
scattering  cross  section  can  be  though  of  as  the  probability  that  a 
scattering  event  will  occur  and  couples  the  backscattered  yield  with 
the  quantitative  analysis.  The  differential  scattering  cross  section 
including  the  cor-ection  for  target  recoil  is: 


do  (6) 
dft 


zz^2  2 (cose  + [i-(^  s,„  e)21  }2 

2ESin  20  [l-(^  SinS)2]*5 


(A.  2) 


The  average  value  of  the  differential  scattering  cross  section,  o(0), 
over  the  solid  angle  subtended  by  the  detector,  ft,  is: 


do  (6) 


(A.  3) 


Exceptions  to  the  Rutherford  formula.  Equation  A. 2,  occur  when  the 
energy  of  the  incident  particle  is  low  and  the  target  is  heavy 
(e.g.  10  KeV  He+  on  Au)  and  when  the  energy  of  the  incident  particle 
is  high  and  the  target  is  light  (in  this  energy  range  reaction 
resonances  occur) . 

The  energy  loss  factor  arises  from  the  fact  that  when  an  energetic 
charged  particle  passes  through  matter  it  will  lose  energy  in  a ca  - 
culable  way  due  to  atomic  ionization  and  excitation  of  the  target 
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atoms.  Figure  A. 4 shows  graphically  how  the  observed  energy  can  be 

associated  with  the  depth  at  which  the  scattering  event  occurred. 

If  an  alpha  particle  is  incident  on  the  sample  surface  with  energy 

E (0)  it  will  have  energy  E (x)  at  some  distance  x beneath  the  sample 
a a 

surface.  At  some  distance,  say  x^,  beneath  the  surface  the  particle 

will  scatter  from  a target  nucleus,  assumed  to  be  type  1,  and  the 

scattered  particle  will  have  energy  K^E^(x)  or  Eo^x^.x^).  The 

scattered  particle  will  also  lose  energy  as  it  traverses  the  sample 

and  will  be  observed  to  have  energy  Ea^(0,x^).  A9  can  be  seen  there 

will  be  a one-to-one  correlation  between  observed  energy  and  the 

scattering  depth  for  each  component. 

Experimental.  A scattering  chamber  utilized  for  RBS  is  shown  in 

4 + 

Figure  A. 5.  Typically,  2.0  MeV  He  are  used  as  incident  particles 
to  minimize  Coulomb  scattering  and  resonant  scattering. 

Applications . The  applications  of  RBS  are  far  reaching  (72)  due 
to  the  ability  to  measure  composition  variations  and  impurity  dis- 
tributions as  a function  of  depth  with  a depth  resolution  of  10  nm 
over  the  first  300  nm  of  sample  without  the  need  for  serial  section- 


ing or  sputtering.  The  other  advantage  of  RBS  is  its  ability  to 


perform  quantitative  depth  analysis  without  the  need  for  standards 


of  similar  composition,  the  technique  is  is  absolute . A few  of  the 


many  applications  of  RBS  will  be  discussed  below 


One  of  the  most  important  applications  of  RBS  has  been  in  the 


area  of  thin  film  studies  (72-79).  Figure  A. 6 (74)  shows  the  power 


of  the  technique.  The  alpha  spectrum  shown  is  for  a multilayer 


j 


F 
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sample  of  ZnS  aiid  ThF^  on  quartz,  there  are  twelve  105.6  nm  layers 
of  ThF^  separated  by  67.0  nm  layers  of  ZnS.  Due  to  the  high  differ- 
ential cross  section  for  Thorium,  the  peaks  associated  with  zinc, 
sulfur,  and  Fluorine  cannot  be  directly  seen.  The  arrows  indicate 
the  position  the  peaks  should  occur  at  if  located  at  the  sample  sur- 
face . 

Another  important  application  of  RBS  has  been  in  the  area  of  thin 
film  stoichimetry  (76-79).  Wenz  and  Hoffman  (78)  and  Wenz  (79)  have 
studied  the  intrinsic  stress  and  microstructure  as  a function  of 
stoichimetry  of  ZnS^  and  TiO^  films.  Figure  A. 7 (79)  shows  a compari- 
son of  a ZnS^  film  deposited  on  a nickel  substrate  (solid  line)  with 
that  of  a ZnSx  film  of  the  same  thickness  deposited  on  a carbon  sub- 
strate (dots). 

Other  applications  of  RBS  to  materials  related  problems  have  been 
in  the  area  of  solid  state  diffusion  (80),  oxidation  mechanism 
studies  (81-86),  the  sputtering  process  (87)  and  surface  roughness 
(88).  Schmid  and  Ryssel  (88)  have  measured  surface  roughness  of  the 
order  of  300  nm  by  observing  the  front  edge  of  the  observed  spectra. 

By  using  resonant  backscattering  Mackenzie  and  Armitage  (89)  have 
measured  sample  perosity.  They  used  the  2.66  MeV  resonance  in  the 
proton  scattering  cross  section  in  oxygen-16  to  measure  pore  sizes 
as  small  as  1.0  nm  by  obersving  the  resonance  width.  The  observed 
resonance  width  will  depend  on  the  number  and  size  of  the  pores  en- 
countered before  and  after  resonant  backscattering  . 
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A.  4 Channeling 

Introduction.  Historically  channeling  was  first  suggested  by 
Stork  and  Wendt  (90)  in  1912  for  single  crystal  studies.  Of  all  the 
nuclear  microanalytical  techniques  it  is  the  most  difficult,  but 
has  the  ability  to  determine  impurity  and  defect  locations  to  an 
accuracy  of  M.0  pm  (91) . 

Concept . Channeling  involves  directing  the  incident  beam  along 
a crystallographic  axis  so  the  interaction  between  the  beam  and 
matrix  is  small.  Scattering  events  will  occur  when  an  impurity  or 
defect  interacts  with  the  channeled  beam,  the  energy  of  the  scattered 
particle  is  then  a measure  of  the  location  of  the  scattering  center. 

Experimental.  The  experimental  setup  for  channeling  experi- 
ments is  similar  to  that  for  RBS,  with  the  exception  that  a better 
collimated  beam  is  required  (<0.05°)  and  two-axis  goniometer  is 
required  for  precise  sample  positioning  (<0.05°).  It  is  also 
necessary  to  be  able  to  move  the  beam  spot  laterally  to  check  for 
radiation  damage  induced  by  the  incident  and  scattered  particles. 

It  may  also  be  desirable  to  position  the  detector  to  allow  both 
channeling  and  blocking  (channeling  of  the  scattered  particle) 
kown  as  double  alinement. 

Applications.  The  most  common  uses  of  channeling  to  date  have 
been  in  the  areas  of  impurity  distributions  (80,92,93),  defect 
studies  (94,95),  and  radiation  damage  (96,97).  An  example  of  the 
type  of  information  that  can  be  obtained  by  channeling  studies  is 
shown  in  Figure  A. 8 (93).  The  spectra  are  from  an  as-grown  silicon 
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A. 8 RBS  spectra  for  1.8  MeV  He  incident  on  a 6 x 10^"*  As-m  doped 
Si  crystal. 
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crystal  containing  6 x 10  As-m  . The  solid  circles  are  for  a 
random  direction,  the  triangles  are  for  a crystal  alined  in  the 
<110>  direction,  and  the  open  circles  are  for  an  alined  sample 

-2 

in  the  <110>  direction  after  a random  bombardment  of  6.5  pC-mm 
of  1.8  MeV  alpha  particles. 

Recently,  channeling  and  blocking  have  been  combined  to  locate 
surface  atoms  with  a precision  of  2 pm  to  study  surface  relaxation 
(98-100). 


A. 5 Particle  Induced  Nuclear  Reactions 

Introduction.  When  it  is  desirable  to  obtain  concentration 
versus  depth  profiles  of  a specific  isotope  it  is  often  possible  to 
use  particle  induced  nuclear  reactions.  There  are  two  possible 
types  of  particle  induced  nuclear  reactions,  namely  1)  the  creation 
of  a radioactive  species  and  c2-  ^ ^ >•  direct 

observation  of  a nuclear  reaction. 

Creation  of  a radioactive  species  involves  the  conversion,  by 
nuclear  reaction,  of  the  isotope  of  interest  to  a radioactive  species 
(101-104) . Depth  profiling  is  then  accomplished  by  sectioning  the 
specimen  and  the  determing  the  concentration  of  radioisotope  by 
the  section  activity.  Because  of  the  similarity  to  the  conventional 
radiotracer  experiments  after  the  sample  has  been  Irradiated,  no 
further  discussion  of  induced  radioactivity  will  be  included. 

The  direct  observation  of  nuclear  reactions  involves  bombard- 
ing a sample  containing  the  Isotope  of  interest  with  a monoenergetic 
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crystal  containing  6 x 10  As-m  . The  solid  circles  are  for  a 
random  direction,  the  triangles  are  for  a crystal  alined  in  the 
<110 > direction,  and  the  open  circles  are  for  an  alined  sample 
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in  the  <1 10 > direction  after  a random  bombardment  of  6.5  pC-run 
of  1.8  MeV  alpha  particles. 

Recently,  channeling  and  blocking  have  been  combined  to  lo- 
cate surface  atoms  with  a precision  of  2 pm  to  study  surface  re- 
laxation (98-100). 

A. 5 Particle  Induced  Nuclear  Reactions 

Introduction.  When  it  is  desirable  to  obtain  concentration 
versus  depth  profiles  of  a specific  isotope  it  is  often  possible  to 
use  particle  induced  nuclear  reactions.  There  are  two  possible  types 
of  particle  induced  nuclear  reactions,  namely  1)  the  creation  of  a 
radioactive  species  and  2)  the  direct  observation  of  a nuclear  reac- 
tion. 

Creation  of  a radioactive  species  involves  the  conversion,  by 
nuclear  reaction,  of  the  isotope  of  interest  to  a radioactive  species 
(101-104).  Depth  profiling  is  then  accomplished  by  sectioning  the 
specimen  and  the  determining  the  concentration  of  radioisotope  by  the 
section  activity.  Because  of  the  similarity  to  the  conventional  radio- 
tracer  experiments  after  the  sample  has  been  irradiated,  no  further  dis- 
cussion of  induced  radioactivity  will  be  included. 

The  direct  observation  of  nuclear  reactions  involves  bombarding  a 
sample  containing  the  isotope  of  interest  with  a monoenergetic 


beam  of  charged  particles  and  observing  the  instantaneously  emitted 
energetic  particles  (2,38,105-143).  The  remainder  of  this  section 
will  be  devoted  to  this  technique. 

General.  In  general  the  direct  observation  of  nuclear  reactions 
has  the  following  properties  (105) ; 

(a)  The  technique  is  restricted  to  the  first  10  pm  of  the  sample. 

(b)  With  most  reactions  there  is  no  natural  background. 

(c)  Most  of  the  reactions  have  a high  Q value  (Q  is  the  energy 
equivalent  to  the  mass  difference  between  the  reactant  and 
product  nuclei)  which  allows  the  use  of  low  energy  particles, 
thus  reducing  local  heating  of  the  target  and  avoiding  com- 
peting reactions. 

(d)  Coulomb  barriers  prevent  the  reaction  of  most  nuclei  larger 
than  chlorine  with  the  low  energy  particles  employed,  thus 
allowing  the  determination  of  light  elements  in  heavy 
matrices . 

(e)  The  nuclear  reactions  employed  are  very  specific,  in  that 
two  isotopes,  even  of  the  same  element,  will  react  quite 
differently. 

(f)  The  results  obtained  are  independent  of  the  matrix  in  which 
the  isotope  is  embedded,  since  the  physical  or  chemical 
state  of  the  atom  does  not  effect  the  nuclear  reaction. 

(g)  The  results  are  quantitative. 
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(h)  The  sensitivity  is  high,  in  that  as  little  as  10  kg-m 
can  be  determined. 

(i)  It  is  possible  to  measure  concentration  versus  depth  pro- 
files. 

A partial  list  of  isotopes  which  can  be  used  for  the  direct 
observation  of  nuclear  reactions  is  given  in  Table  A. 2 along  with 
the  Q value  of  the  reaction. 

Theory . When  a monenergetic  beam  of  charged  particles  of 
energy  E^O)  is  incident  on  a sample  the  particles  will  lose  energy 
almost  linearly  as  they  traverse  the  sample  and  will  have  energy 
E^(x)  at  any  position  within  the  sample,  see  Figure  A. 9.  At  some 
depth,  say  x^»  beneath  the  surface  the  particle  will  react  with  a 
nucleus  of  the  isotope  of  interest  and  will  emit  a particle  of 
energy  Ee(x^,x^).  The  energy  of  the  emitted  particle  can  be  cal- 
culated from  the  following  equation  (145) 


E (x,x) 
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(m.m  E (x)) 

r * 1 1 


m + m 
e r 
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' e r 
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m + m 
e r 


(A. 4) 


Where  m^,  m^,  and  m^  are  the  masses  of  the  incident  particle,  the 
emitted  particle,  and  the  residual  nucleus  respectively.  The  emitted 
particle  will  also  lose  energy  almost  linearly  as  it  traverses  the 


TABLE  A . 1 


Q Values  for  Some  Nuclear  Reactions  Induced  by  Charged 
Particles . 
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sample  and  will  have  energy  Ee(0,x^)at  the  sample  surface.  As  can 
be  seen  from  Figure  A. 9 there  will  be  a one-to-one  correspondence 
between  the  observed  energy  of  the  emitted  particle  and  the  depth  at 
which  it  was  created. 

The  number  of  particles  detected  in  an  energy  interval,  N£)(xi), 
associated  with  mean  depth  x.^  is  given  by  the  following  equation: 

NQ(xi)  «=  Nx(x^) jo (x^)AXjAfl  (A. 5) 


where  NC(x^)  is  the  number  concenteration  of  isotope  with  mean  depth 
x^,  j is  the  integrated  flux  of  particles,  ©(x^  is  the  differential 
cross  section  at  mean  depth  x , Ax^  is  the  depth  interval  and  An 
is  the  solid  angle  of  observation.  It  can  therefore  be  seen  that 
ideally  the  number  of  counts  recorded  in  a given  depth  interval  will 
be  a measure  of  the  concentration  of  isotope  in  the  interval. 

In  reality  the  ideal  spectrum  will  be  convoluted  by  a spreading 
function  P(E^,Ej),  which  is  the  probability  that  a particle  that 
should  be  recorded  in  energy  interval  i will  be  recorded  in  interval  j. 
The  distorted  spectrum,  [Nq]  is  then  given  by: 


Experimental . Experimentally  a sample  is  irradiated  with  a mono 
energetic  beam  of  particles  and  the  energetic  particle  emitted  at  a 
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fixed  laboratory  scattering  angle  observed  using  a conventional 
solid  state  detector.  Because  the  cross  section  for  Rutherford 


scattering  is  greater  than  ('v-lO  times)  the  reaction  cross  section 
it  is  necessary  to  separate  the  scattered  incident  particles  from 
the  particle  emitted  by  the  reaction  to  prevent  pulse  pile-up  in 
the  electronics.  Two  methods  exist  to  separate  the  scattered  par- 
ticles from  the  emitted  particles,  they  are:  1)  the  placing  of  an 
absorber  in  front  of  the  detector,  typically  mylar  (106-108)  or 
nickel  (109)  and  2)  the  use  of  a magnetic  spectrometer  (33,110)  to 
separate  the  particles  according  to  their  momentum.  Figure  A. 10 
(33)  shows  the  advantage  of  using  the  magnetic  spectrometer  over 

absorbers.  Figure  A. 10. a is  a typical  spectrum  using  a mylar  ab- 

18 

sorber,  the  sample  was  polycrystalline  A^O^  annealed  in  21%  0£ 

4 

for  8.64  x 10  seconds  at  1641  K.  The  reaction  induced  is  the 
18  15 

0(p,a)  N reaction  with  an  incident  proton  energy  of  0.75  MeV. 
Figure  A.IO.b  is  of  the  same  sample  only  using  the  magnetic  spec- 
trometer, all  other  parameters  are  the  same.  The  difference  in 
alpha  yield  is  due  to  a difference  in  solid  angle. 

Figure  A. 11  shows  a typical  scattering  chamber  for  the  direct 
observation  of  nuclear  reactions. 

Analysis  of  the  experimental  data  can  be  performed  by  either 
the  resonance  technique,  where  the  concentration  profile  is  obtained 
by  walking  an  isolated  resonance  through  a sample  as  was  done  in 
this  discertation  or  by  the  single  spectrum  technique  where  the 
entire  spectrum  is  deconvoluted  to  obtain  the  concentration  versus 

depth  profile  first  used  by  Robin  (2). 


tab. 
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reaction  (33). 


A. 11  Typical  scattering  chamber  for  the  direct  observation  of  nuclear  reactions. 


Applications . The  direct  observation  of  nuclear  reactions 
has  three  general  areas  of  application  namely:  1)  trace  concentra- 
tions, 2)  diffusion,  aid  3)  ccrrosion.  Trace  concentrations  of 

UB  (120),  12C  (121),  UN  (122-124) , 160  (124-126),  180  (120,127,128), 
23 

and  Na  (130)  have  been  detected  as  a function  of  depth  by  the 

direct  observation  of  nuclear  reactions.  North  and  Lightowlers  (120) 

1L  8 18  15 

utilized  the  TKp.a^)  Be  -+  and  0(p,a)  N reactions  to 

study  surface  contamination  due  to  solvents.  Figure  A. 12  shows  the 


alpha  spectra  obtained  from  a control  sample  (solid  circles)  and 

from  a sample  cleaned  in  electronic  grade  methanol  (open  circles). 

Average  concentrations  for  the  control  sample  are  (2±l)x  10^8 
-2  19  -2 

B-m  and  (1.4±0.3)  x 10  0-m  . The  corresponding  concentrations 

16  _2 

for  the  methanol  contaminated  sample  are  (4.0±0.4)  x 10  B-m  and 
(5. 1+0. 7)  x 1019  0-m~2. 

The  second  area  of  application  has  been  in  the  field  of  solid 

state  diffusion  processes  (2,  106-108,  131-137).  Figure  A. 13  (136) 

18 

shows  the  variation  in  the  observed  alpha  spectra  from  the  0(p,a) 
^"*N  reaction  as  the  value  of  Dt  is  varied.  The  incident  proton 
energy  is  0.75  MeV. 

The  third  area  where  the  direct  observation  of  nuclear  reactions 
has  proven  useful  is  in  oxidation  studies  (38,110,114,138,139). 

Figure  A.  14  (38)  shows  the  type  of  data  that  can  be  obtained  from 


oxidation  studies.  The  alpha  spectrum  is  from  a silicon  coupon 
which  was  thermally  oxidized  in  ^^2  at  *295  K to  produce  an  0.5  pm 
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esence  of  two  process 


TABLE  A. 3 


Possible  Coincidence  Measurements 
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A. 6 Current  Trends 


The  trends  in  nuclear  microanalysis  are  to  enhance  both  the 

depth  and  lateral  resolution  for  surface  analysis.  Williams  (39) 

has  suggested  using  low  angle  incident  beams  to  improve  the  depth 

resolution  in  RBS  and  particle  induced  nuclear  reactions.  By  using 

an  incident  angle  of  5°  from  the  sample  surface  and  using  a solid 

state  detector  (15  keV  FWHM)  positioned  at  168°  it  is  possible  to 

obtain  a depth  resolution  of  2-3  nra  in  RBS  experiments. 

Pierce  and  Duddleston  (145)  have  reduced  the  beam  size  to 
2 

25  pm  and  used  a two  dimensional  scanning  technique  to  map  light 

elements.  By  combining  PIXE,  for  iron  determination,  and  the  direct 
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observation  of  nuclear  reactions,  C(d,p)  C,  0(d,p)  0,  and 

^B(d,p)^*B,  they  were  able  to  map  the  concentration  distributions 
of  iron,  carbon,  oxygen,  and  beron  across  a weld. 

A. 7 Comparison  of  Techniques 

Recently  several  authors  (146-149)  have  compared  the  results 
obtained  by  nuclear  techniques  with  those  obtained  by  Auger  Electron 
Spectroscopy  (AES),  and  X-Ray  Photoelectron  Spectroscopy  (XPS) . 

Table  A. 4 gives  a comparison  of  the  various  techniques  (146,149). 

One  should  bear  in  mind  that  none  of  the  techniques  are  them- 
selves the  answer  to  all  surface  analysis  problems,  but  with  a 
combination  of  techniques  it  is  possible  to  obtain  a complete  chemi- 
cal analysis. 
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APPENDIX  B 

THE  15N(p,ari)12C  CROSS  SECTION 


B.l  Introduction 


When  nitrogen-15  is  bombarded  with  low  energy  protons,  <2  MeV, 
it  is  energetically  possible  for  the  following  three  reactions  to 


15N  + p 12C  + a 

o 

15n  + p , 12c*  + .J 

15N  + p 16o  + Y 


These  reactions  have  been  studied  by  many  investigators  (150-159) 
with  one  of  two  objectives.  The  first  of  these  was  to  gain  insight 
into  the  nuclear  energy  levels  of  the  ^0  compound  nucleus  and  the 
second  was  to  study  the  carbon-nitrogen  cycle. 

When  Bethe  (160,161)  investigated  the  carbon-nitrogen  cycle 
as  a source  of  stellar  energy,  the  cross  section  for  these  reactions 
became  of  interest  at  stellar  energies,  near  30  keV.  The  burning 
of  hydrogen  in  second  generation  stars  with  masses  greater  than  that 
of  the  sum  in  predominantly  through  the  carbon-nitrogen  cycle  which 
involves  the  first  two  reactions  given  by  (159)  : 

12C(p,y)13N(B+v)13C(p,Y)14N(p,y)150(B+v)15N(p,a)12C 


*?—  - .i  _ 


This  sequence  of  reactions  releases  approximately  25  MeV  of  energy 
per  cycle.  The  third  reaction  allows  for  leakage  from  the  cycle, 
however,  it  is  returned  to  the  cycle  by  one  of  the  following 
sequences  (159): 


or 


15N(p,y)160(p,y)17F(B  v)170(p,a)1<!,N 


15N(p,Y)160(p,Y)17F(e+v)170(p.Y)18F(B+v)180(p,a)l5N 


It  is  reported  by  Rolfs  and  Rodney  (159)  that  the  rate  of  the 

8N(p,y)  0 reaction  will  determine  the  overall  abundance  of  ^0, 

17  18 

0,  and  0 synthesized  in  the  cycle  and  will  therefore  play  an 

important  role  in  stellar  nucleosynthesis. 

This  study  orginated  for  two  reasons.  First,  the  differential 
15  12 

cross  section  for  the  N(p,aQ)  C reaction  was  not  reported  for  a 
laboratory  scattering  angle  of  165°  and  secondly,  due  to  the  author's 
interest  in  nuclear  physics.  Realizing  that  roost  of  the  readers  are 
not  familiar  with  nuclear  physics  a glossary  of  the  symbols  used  in 
this  Appendix  is  given  in  the  next  section. 


B.2  Glossary  of  Symbols 
A = nucleon  number 

E ^energy  of  the  incident  partcle  inappropriate  reference 
frame 
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Eq=  resonant  energy  in  laboratory  reference  frame  (this  term  is  not 
a constant)  defined  as 


E = E.  + A, 
o X X 


(B.l) 


E^  = channel  energy  in  the  laboratory  reference  frame. 

Cn  »P ) = regular  solution  to  the  radial  wave  equation  (regular 
coulomb  function. 

G^q.p)  = irregular  solution  to  the  readial  wave  equation  (irregular 
coulomb  function) . 

1,1'  = intrinsic  spin  of  the  target  nucleus  and  recoil  nucleus  re- 
spectively . 

J = total  angular  momentum  of  the  system  give  by 


J = s + fc  ■=  s'  +«.' 


(s  is  the  channel  spin  defined  below) 


(B . 2) 


JQ  = total  angular  momentum  of  a particular  nuclear  energy  level. 


N = number  density. 


Nq  = number  of  observed  alphas 

PL(Cos  0)  = Legendre  Polynomial  given  by 


PL(Cos  0) 


(St„9)ZL 

2 L!  d(Cos9) 


(B.3) 


Q = the  energy  equivalent  to  the  rest  mass  difference  between  the 
reactant  and  product  nuclei. 
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S = the  scattering  matrix  whose  elements  are  amplitudes  of  prob- 


ability flux  (S  is  both  unitary  and  symmetric) . 

Z = atomic  number. 

g^  = a real  quantity  which  is  related  to  the  partial  width  of  chan- 
nel a,  T , by  : 
a 

ga  ■=  ± (ra)*5  (B.4) 

i,i'  = the  intrinsic  spin  of  the  incident  and  emerging  particles 
respictively . 

j = integrated  proton  flux. 

k = wave  number. 

m^  = atomic  mass  of  species  i. 

r = interparticle  separation. 

r = 1.3  x 10  ^ m. 
o 

r H the  screening  radii  for  the  Coulomb  field  in  channel  a . 
a 

s,s'  = the  incoming  and  outgoing  channel  spins  respectively  given 
by  : 

s = i + I (B.5) 

s - t'  + t'  (B .6) 

i.,1'  = the  orbital  angular  momentum  of  the  incoming  and  outgoing 
channels  respectively. 

v = relative  velocity  of  the  colliding  particles. 
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r = the  total  resonance  width  in  the  laboratory  reference  frame 
which  is  a measure  of  the  mean  lifetime  of  the  nuclear  level. 

A^  = the  level  shift  which  relates  the  channel  energy,  E^  , and  the 
resonant  energy,  Eq,  by; 


E “ E + A. 
o \ X 


(B  .7) 


it  = parity  which  is  a measure  of  symmetry  under  space  inversion. 
ttq=  parity  of  a particular  nuclear  energy  level, 
ft  = solid  angle. 

a,a'=  channel  index  for  the  incoming  and  outgoing  channels  respec- 
tively, and  defines  the  type  and  state  of  the  particles. 

B = the  resonant  phase  shift  which  is  an  angular  measure  of  the 
deviation  of  the  energy  from  the  resonance  energy,  defined  as 


-1  2(E-Eo) 

B - Tan  


(B  .8) 


a term  given  as 


0.1574  Z Z m ■ 
a x a 


(B  .9) 


for  the  incident  channel  and 


n - 0.1574  ZbZy(^),S 


(B.10) 


for  the  exit  channel 


HiMHMBaflikr 


6 = scattering  angle. 


r 


X r deBroglie  wavelength  (divided  by  2n) . 


V reduced  mass. 


= phaseshift  for  potential  (hard  sphere)  scattering  (quantum 


mechanical  scattering  from  a hard  sphere). 


p = parameter  defined  as: 


p = Kr 


(B.ll) 


o = cross  section,  a measure  of  the  probability  of  a reaction 


occurring . 


a = phaseshift  for  Coulomb  scattering  from  an  inpenetrable  sphere. 


B.3  Theoretical  Treatment  of  the  Reaction  Cross  Section 


The  theoretical  treatment  presented  in  this  section  follows 


that  given  by  Blatt  and  Biedenharn  (162)  for  the  reaction: 


a + X -►  Y + b 


(B.12) 


Here  "a"  represents  the  incident  particle  with  intrinsic  spin  i 


which  collides  with  the  target  nucleus  X,  which  has  intrinsic  spin 


I.  Particle  b with  intrinsic  spin  i'  emerges  at  an  angle  0 to  the 


incident  beam  in  the  center-of-mass  system  and  the  residual  nucleus, 


Y recoils  in  the  opposite  direction.  The  intrinsic  spin  of  Y is 


V. 


It  is  convenient  to  describe  the  system  before  and  after  the 


reaction  with  three  parameters:  a(  or  a’),  the  channel  index,  s 


(or  s’),  the  channel  spin,  and  £(or  £’),  the  orbital  angular  momen- 


tum in  the  center-of-mass  system.  The  total  angular  momentum  J,  of 


the  system  is  given  by  vector  addition  of  s and  £ or  s'  and  £' 


S + £ ■=  J « 8'  + £' 


(B .13) 


For  the  case  where  the  reaction  proceeds  via  a discrete  energy  level 
in  a compound  nucleus,  the  total  angular  momentum  of  the  system  is 


equal  to  the  angular  momentum  of  the  energy  level,  J£ 


s + £ = J = s’  + £' 
o 


(B.14) 


The  total  angular  momentum  of  the  system  is  always  conserved  (s 


and  £ are  not  conserved  quantities). 


For  strong  interactions,  which  nuclear  reactions  are,  parity. 


it  , is  also  a conserved  quantity.  If  the  parity  of  the  compound 


nucleus  is  given  by  tto  the  following  condition  must  be  satisfied: 


(-)  It  = IT 

a o 


<->  V 


(B.15) 


where  *a  and  ttq  are  the  channel  parities  for  channel  a and  a' 


respectively  given  by: 


(B.16) 


■ * rat,.  > ximmin-t 


In  general,  one  usually  introduces  a unitary,  symmetric  matrix, 

S,  known  as  the  scattering  matrix  whose  elements  are  amplitudes  of 
probability  flux  for  a collision  where  channel  as l goes  to  channel 
a's'i’  with  total  angular  momentum  J (S  is  a unitary  matrix  in 
order  to  conserve  the  normaliz  tion  of  the  wave  function  with  time 
and  is  symmetric  in  order  to  conserve  time  reversal) . The  cross  section 
is  now  obtained  by  averaging  over  the  incident  spin  directions  and 
summing  over  the  final  spin  directions,  and  is  given  by: 


x K(J1d|«.1;J2J.2£2;s,s;e)dfi  (B.18) 

where  * is  the  deBroglie  wavelength  (divided  by  2n)  and  K is  pure- 
ly geometrical  in  that  it  is  independent  of  the  nature  of  the 
channels  and  the  dynamics  of  the  collision  and  can  be  written  in 
terms  of:  Z which  are  defined  as: 


I ( 

I : 


r 


L-*l+* 2 

Z(11J1I12J2,SL)  " 1 [(241+l)(2l2+l)(2J  +1) 


2(J2+l)]35WU1J1£2J2,sL)ait200|jl1il2L0) 


(B.19) 


is  the  Racah  coefficient  and  | is  the 

Clebsch-Gordon  coefficient  (vector  addition).  The  first  term, 
L‘£l+£2 

i , is  always  real*  (±1) . 

It  is  now  possible  to  write  the  differential  cross  section  as 
a sum  of  terms  in  the  following  manner: 


2 -1 

“ * (2s+l)  l B.  (a's'  ;as)PT  (Cos0)dft 

Cl  b *CXo  (X  _ . L L 

L=0 


where: 


B^(a's' ;as) 


I l l l I,  I, 

J1  J2  £1  £2  £i  £2 


Z(£1J1£2J2’SL)Z(£IJ1£2J2’S’L) 


R.P.[(6  , 5 ' 6o'o  - 6 , *1) 

a as  s a s £ ;asl 


'o^s * s^ £.2R 2 's ' i2;asil2  ^ 


(B.20) 


(B.21) 


I 


The  notation  R.P.[  ] stands  for  the  real  part  of  the  term  in 

brackets.  The  summations  in  Equation  B.21  all  run  from  zero  to 
infinity,  however  only  one  of  the  sums  is  actually  infinite  due  to 
the  selection  rules  for  nonvanishing  Z coefficients.  By  conserva- 
tion of  time  reversal: 


B (a's';as)  = B (asja's’) 

L L 


(B .22) 


Equation  B.21  should  be  contracted  so  that  each  term  is  counted 
only  once.  As  written,  identical  contributions  are  obtained  when 
the  terms  are  interchanged.  The  generalized  and  contracted  form  of 
Equation  B.21  is  given  by: 

s ' -s  «o  J+s  J+s ' 

B.  (a's*  ;as)  - •L~ ’'-r-  III  Z(!J£J.sL) 

^ 4 J=0  £-|J-s|  l^J-s'l 


za’Jt’J.s'L)  VaVsV!  - Sa,s,f.'  ;as£ 


1=1 


s -s 


J^+S 


J^+s ' 


Jj+s  J2+s' 


l . I 


l <_  I ..  . X ...  X 


+ 2 J -0  £ -1^-*!  A|-|J1-s,|  J2-Jx+1  | j2“s I a2=  I J2_S*  ^ 


Z(Vlt2J2,Sl)Z(liJll2J2’S'L) 


R.P.  [ ] 
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. 


Jj+S  Jx+8 ' 

+ I [ Z(UJJ  sL)2(i:J1t;J1,s'L)R.P.[J,«J1] 

£2m<ll+1  42"lJrS’l  1121  1121  2 1 


V 


+ l 


*2"Ari 


Z(VlJllJl’8L)ZaiJlll2Jl,S,L)RPfJ2MJl’£2*)llJ) 


(B.23) 


The  term  R.P.[  ] stands  for  the  corresponding  expression  in 

Equation  B.21.  The  additional  restriction  that: 

+ l2  - L - even;  S,|  + - L = even  (B.24) 

will  also  reduce  the  number  of  terms  in  Equation  B.23. 

In  obtaining  Equation  B.23  no  assumptions  have  been  made  about 
the  reaction  mechanism,  with  the  exception  that  two  particles  enter 
and  two  particles  emerge.  For  any  specific  case  many  simplifications 
can  be  made,  such  as,  the  case  where  the  reaction  proceeds  via  a 
definite  resonance  level  of  a compound  nucleus  with  angular  momen- 
tum J and  parity  ir  . Blatt  and  Biedenham  (162)  now  find  it  con- 
o o 

venient  to  introduce  a real  quantity,  g ^ » which  is  related  to 
the  partial  width,  I*  g,  by  the  following: 


The  ambiguity  In  sign  does  not  exist  In  the  formula  for  the  cross 

section,  however  a study  of  the  angular  distribution  for  the  reaction 

allows  the  determination  of  both  the  magnitude  of  T and  the  sign 

as  x 

of  g . . 
eas  l 

It  Is  also  necessary  to  introduce  the  phaseshlfts  for  potential 
Chard  sphere)  scattering,  £ , which  is  given  by  Equation  B.26  for 

charged  particles: 


exp(2Ut) 


exp(2iafc) 


G^Cn.Pi-iF^Cn.p) 
G^Cn.p)  +iFt(n,p) 


(B.26) 


where  F^(n,p)  and  G^(n,p)  are  the  regular  and  irregular  solutions  to  the 
radial  wave  equation  outside  the  nuclear  surface.  For  the  incident 
channel  the  value  of  n is  given  by: 


n 


0.1574  Z Zvm  ■ 
a X a 


(B.27) 


where  Zfl  and  Z^  are  the  atomic  numbers  of  the  incident  particle 
target  nucleus  respectively,  m is  the  mass  of  the  incident  parti- 
cle  in  amu  and  E is  the  incident  particles  energy,  in  MeV,  in  the 
laboratory  reference  frame.  For  the  exit  channel  n is  given  by 
(163): 


n 


0.1574  ZfeZY  r )?5 


(B.28) 


where 


and 


are  the  atomic  numbers  of  the 
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emitted  particle  and  recoil  nucleus  respectively,  p'  is  the  reduced 
mass  given  by: 


Vy 


and  E'  is  the  energy  available  to  the  reaction  given  by: 


v 

E’  - ; + Q 


The  parameter  p is  given  by: 


p = kr 


(B .29) 


(B . 30) 


(B.31) 


where  u is  the  reduced  mass,  v is  the  relative  velocity  of  the  col- 
liding particles,  and  r is  the  interparticle  separation.  It  is  also 

possible  to  approximate  p by: 

1/3  wo 


1/3 

ro(Aa  + AX  > 


(B.32) 


where  rQ  is  approximately  1.3  x 10  m and  A is  the  nucleon  number. 

The  phaseshift  for  Coulomb  scattering  from  an  impentrable 
sphere,  , is  given  by  the  following  formula: 


*>xd  f 2io  ) ■=  Odin)  1*0 )_L 

expuicy  (t-in)(i_i-in). ..  (i-in)  (-in)! 


(B.33) 


i 


jrr. 


The  observed  resonant  energy,  Eq,  is  not  a constant,  but  is 
dependent,  on  the  channel  energy  through  a quantity  . The  resonant 
energy  should  be  given  by: 

E - + A,  (B . 34) 

O X A 

However  in  practice  E^  can  be  considered  a constant. 

The  scattering  matrix  elements  vanish  unless  J = J and  it  *=  r 

o o 

for  the  reaction  cross  section.  The  scattering  matrix  can  now  be 
written  in  the  form: 


S , ,.i  .‘Vo  * exp[i(£  - n ln2k  r )]exp[i(£  Iol-n  ,ln2k  ,r  ,)] 
o's  1 ;asl  K oi.  a a a r sa  Jt'  a a a 


X ‘VaVaV*  +i  <for  J " V * = V 

o 


(B.35) 


Ilere  ra  and  r , are  the  screening  radii  for  the  Coulomb  field  it. 
channels  a and  a'  respectively.  In  the  final  equation  the  screen- 
ing radii  will  drop  out.  The  total  width  for  the  level  J , it  , E 
° o o o 

is  given  by: 

r * l l l r , (B.36) 

L L „ as£ 
a s f. 

and  depends  on  the  energy.  Equation  B.35  assumes  that  resonances  of 
different  J and  it  are  far  enough  away  in  energy  that  their  influence 
can  be  neglected. 
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Substituting  Equation  B.35  into  Equation  B.21  and  allowing: 


BL(a's';as)  «=  ^(a's'jas) 


(B .37) 


gives : 


J +s  J +s  J +s ' 

o o o 


(a 's  ' ;as) 


4[(E-Eo)2+(Hr)2]  fi-1=|ao-s|  i2*|Jo-s|  i;-|Jo-s'|  t'-|Jo-s’| 


J +s* 
o 

l 


* Z(tlJot2Jo,sL)ZaiV2Jo,S'l)  8as£1  8as*2  Vs’^  Vs’*’ 


cosI««£1  ' + 


(for  a»s  + a’ ,s’) 


(B .38) 


The  condition  of  parity  conservation.  Equation  B.15  must  also  be 
imposed  on  Equation  B.38.  As  was  the  case  with  Equation  B.21  each 
term  in  Equation  B.38  will  occur  twice,  it  therefore  should  be  con- 
tracted along  the  same  lines  as  used  to  obtain  Equation  B.23. 

The  total  cross  section,  , can  be  written  as: 


4tt* 

°o(a  'S ' *as)  -2^1  “o 


B (a  ' s ' ; as ) 


it* 


CO  J+S  J+S  ' 

2s+l  J=0  £=|j-s|  i'-lLs'l  (2J+1) 


^a'a^s  's^i.  'l  's ' ;as*. 


(B.39) 


or  substituting 


a (a's'.as) 
o 


9 2J  +1 
ttX2  V5—  I l 

a 2s+l  J {- 


r r 

ast  a's't1 

(e-e  ) 2 + (hr)2 

o 


(for  a,s  j a's') 


(B .40) 


The  sums  in  Equation  B.40  are  the  same  as  those  in  Equation  B.38. 

As  can  be  seen  a measurement  of  the  total  cross  section  does  not 
allow  a determination  of  the  signs  of  the  g's.  However,  in  principle 
a measurement  of  the  differential  cross  section  will  allow  the  signs 
of  the  g's  to  be  determined. 

The  selection  rules  for  nonvanishing  Z coefficients  force 
B to  vanish  if  L violates  one  of  the  following: 


L < 2 *J 


(8.41) 


L * even 


(B.42) 


i;T 


L < 2 l 


max 


(B .43) 
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L < 2 l' 


(B  .44) 


The  formula  given  to  this  point  have  been  for  reactons  which 
have  definite  channel  spins,  however,  the  channel  spins  are  not 
measured  and  the  observed  cross  section  is  fora-*-  a'.  The  differential 
cross  section  is  given  by: 


dV;„  ■ j0  »L<‘‘,;°Vc°s6,i,r! 


(B.45) 


where 


B^Ca'  ;a)  =1^(0'  ;a) 


I+i 

s-ll-il 


l'+i' 

I R.  (a's'  ;as) 

. I T >_■!  • I X 


(B .46) 


It  should  be  emphasized  again  that  Equation  B.45  applies  to  reactions 
only,  that  is  a ^ a* , and  assumes  the  resonances  of  different  J and  w 
are  far  enough  away  in  energy  that  their  influence  can  be  neglected. 

In  practice  the  differential  cross  section  at  an  energy  inter- 
mediate to  two  isolated  resonances  is  given  by: 


do  “ dOj  + do^  + do  £ 


(B .47) 


w.iere  do^  and  do2  are  the  cross  sections  given  by  Equation  B.45  for 
resonances  1 and  2 respectively  and  do^2  is  the  interference  between 
the  two  resonances.  The  interference  term  is  given  by: 


I 


max 


d°12  > * (2i+l)  (21+1)  {4[(E-E  ): 

a , a LSU  1 2 


2,,-% 


+ ftr2>  » Sa.„.  8a'stJ  Z<V<»1t2J02-,U 


Z(tiJol1iJo26'I-)C“  «altl-  «<,*  + ^ + 8!  ' 82> 


PL(Cos  6) 


(B.48) 


The  subscripts  1 and  2 in  the  above  equation  refer  to  resonances  1 
and  2 respectively  and  8 is  an  angular  parameter  which  measures  the 
deviation  of  the  energy  from  the  resonance  energy.  The  value  of  8 
is  given  by: 

, 2(E-E  ) 

8 = Tan  p2-  (B.49) 


The  summation  over  L in  Equation  B.48  runs  from  0 to  L and  includes 

max 


odd  values  of  L if  ir^  f ir^. 
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B.4  Calculation  of  the  N(p,ao)  C Differential  Cross  Section 

15  12 

The  differential  cross  section  for  the  N(p,a  ) C reaction  be- 

o 

tween  0 and  1.5  MeV  can  be  written  as: 
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datotal“  d°338  + d°338/1028+  d°1028  + d°1028/1050  + d°1050  + 


d° 1050/1210  + d°1210 


(B.50) 


where  do^  is  the  differential  cross  section  do  to  the  resonance  at 

energy  i given  by  Equation  B.45  and  do^^  is  the  interference  term 

due  to  resonances  at  energies  i and  j given  by  Equation  B.49.  Each 

term  was  calculated  using  the  resonance  parameters  given  by 

Aj zenberg-Selove  (164).  These  parameters  have  been  tabulated  in 

Table  B.l.  The  partial  width  for  the  a channel,  T , for  the 

o Ot 

O 

1210  keV  resonance  was  calculated  using  Equation  B.36  and  found  to 

be  10.2  keV  (22.5-4.1-8.2  - T ).  Ad'yasevich,  et  al.  (157)  have 

a°  2 

reported  that  the  product  T r to  be  500  keV  for  the  1050  teV 

P “o 

resonance  with  the  terms  r »F  , r,  , equal  to  either  1,500,501  or 

P a lab 
o 

2,250,252  in  keV  respectively.  Calculations  were  carried  out  using 
both  sets  of  parameters. 

In  order  to  determine  the  orbital  angular  momentum  for  each 
channel  it  is  necessary  to  know  the  intrinsic  spins  and  parities  for 
each  of  the  particles  Involved.  For  the  present  case  in  * 

I11  = % , i'11"  0+,  and  I,1T=  0+.  Using  Equations  B.5  and  B.6  the 
channel  spins  are  found  to  be  s = 0 or  1 and  s'=  0.  If  the  1210 

it 

keV  resoannce  is  used  as  an  example  ° = 3 and  to  satisfy  the 
conditions  given  in  Equations  B.14  and  B.15  we  find  that  for  the  in- 
coming channel  s *=  1 and  I = 2 or  4 (s  / 0 by  parity  considerations) 
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TABLE  B.l 
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Resonance  Parameters  for  the  N(p,a0)  C Reaction 


E (keV) 

P 

r (keV) 

P 

r (keV) 

a 

0 

r (keV) 
lab 

J 

0 

338 

1.2 

95 

96 

r 

1028  ± 10 

100 

40 

140  ± 10 

T 

1050  ± 150 

r„r  = 
P “o 

500  keV2 

2+ 

1210  ± 3 

4.1 

resonant 

22.5  ± 1 

3" 
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and  for  the  outgoing  channel  s'  = 0 and  £*  *=  3.  Penetration  con- 
siderations indicate  that  l ■ 4 will  at  most  be  a small  correction 
to  the  cross  section  at  these  energies  and  can  therefore  be  ignored. 

The  Z coefficients  used  in  the  calculations  were  obtained  from 
Biedenharn's  tabulation  of  Racah  Coefficients  (165). 

Figures  B.l  - B.7  show  the  contributions  of  each  of  the  terms  in 
Equation  B.50  at  a laboratory  scattering  angle  of  165°.  It  should 
be  noted  that  because  the  signs  of  the  g's  (±(r)2)  are  not  known  and  the 
interference  terms  may  add  as  the reciprical  of  the  values  given.  After 
each  of  the  individual  terms  were  obtained  they  were  combined  ac- 
cording to  Equation  B.50.  The  possible  total  differential  cross 
sections  at  a laboratory  scattering  angle  of  165°  are  shown  in 
Figures  B.8-B.15.  The  curves  shown  in  parts  a and  b are  for  T ^ 
of  the  1050  keV  resonance  equal  to  252  keV  and  501  keV  respectively. 

The  plus  sign  indicates  that  the  term  added  as  shown  in  Figures  B.l- 


B.7. 
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B.5  Measurement  of  the  NCp.a^l  C Differential  Cross  Section 

15  12 

In  order  to  experimentally  measure  the  N(p,aQ)  C differential 

cross  section  at  a laboratory  angle  of  165°  three  Si^N^  targets 

were  prepared  by  C.  J.  Mogab  (Bell  Laboratories,  Murray  Hill,  New 

Jersey).  The  target  thicknesses  were  0.33  keV(6.5  nm) , 1.30  keV 

(26.0  nm) , and  1.95  keV  (39.0  nm)  to  1.0  MeV  protons  and  the  nitrogen- 

28  _3 

15  density  was  5.5  x 10  atoms-m 


i.36 


Contribution  to  the  differential  cross  section  from  the  338  keV  resonance  (do(338)) 


Contribution  to  the  differential  cross  section  from  the  interference  between  the  338  and  1028 
keV  resonances  (do (338/1028) ) . 


200 


in 
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B.3  Contribution  to  the  differential  cross  section  from  the  1028  keV  resonance  (do(1028)). 


Contribution  to  the  differential  cross  section  from  the  interference  between  the  1028  and  1050 
(T-  252)  keV  resonances  (d  (1028/1050a) ) . 


,5a  Contribution  to  the  differential  cross  section  from  the  1050  (r  - 252)keV  resonance 
(do (1050a) ) . 


B.6a  Contribution  to  the  differential  cross  section  from  the  interference  between  the  1050 
(T  » 252)  and  1210  keV  resonances  (do (1050a/1210) ) . 
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B.6b  Contribution  to  tne  differential  cross  section  from  the  interference  between  the 
1050  (T  ■ 501)  and  1210  keV  resonances  (do (1050b/1210) ) . 


ial  cross  section  from  the 


Possible  differential  cross  section  given  by  do(338)  + do(338/1028) 

+ da (1028)  + da (1028/1050b)  + da(1050b)  + do  (1050b/1210)  + do(1210). 


B.9a  Possible  differential  cross  section  given  by  do(338)  - do(338/1028)  + do(1028) 
+ do (1028/1050a)  + do(1050a)  + do (1050a/1210)  + do(1210). 


+ do (1028/1050b)  + do  (1050b)  + do  (1050b/1210)  + do  (1210). 


ssible  differential  cross  section  given  by  do(338)  + do(338/1028)  + do(1028) 


da (1028/1050b)  + da(1050b)  + da (1050b/1210)  + da  (12] 


•(OUT)  DP  + (OT3!  1 COT)op  - (EOSOT)op  + (B0S0T/830I)  °P 


da(l028/1050a)  + da (1050a)  - da (1050a/1210)  + do (1210). 


B.12b  Possible  differential  cross  section  given  by  do(338)  - do(338/1028)  + do(1028) 


da(1028/1050a)  + da  (1050a)  + da (1050a/1210)  + do (1210). 


da( 1028/1050b)  + do(1050b)  + do (1050b/1210)  + do(1210). 


do (1028/1050a)  + do(1050a)  - do (1050a/1210)  + do(1210). 


do (1028/1050b)  + do(1050b)  - do (1050b/1210)  + do(1210). 


do(1028/1050a)  + do(1050a)  - da (1050a/1210)  + do(1210). 


.15b  Possible  differential  cross  section  given  by  do(338)  - do(1338/1028)  + do(1028) 


lonoenergetic  proton  beams  (±1  keV) , between  0.7  MeV  and  1.5  MeV 
supplied  by  a 3.0  MeV  Van  de  Graaff  accelerator  were  used  to  obtain 
yield  curves  for  each  of  the  three  samples.  The  typical  beam  current 
was  0.4  pA  and  the  beam  profile  was  elliptical  with  a major  axis  of 
3 mm  and  a minor  axis  of  1 mm.  The  energetic  alpha  particles  were 
detected  at  a laboratory  scattering  angle  of  165°  with  a partially 
depleted  surface  barrier  detector  (25  keV  o resolution,  450  mm 
active  area,  and  100  pm  depleted  depth)  which  subtended  a solid 
angle  of  7.24  x 10  ^ sr. 

The  yield  curves  were  used  to  calculate  the  differential  cross 
section  using  the  following  relationship: 


(B.51) 


-where  Nq  is  the  alpha  yield,  N is  the  nitrogen-15  density,  j is  the 
integrated  proton  flux,  x is  the  sample  thickness,  and  Afi  is  the  solid 
angle  of  observation.  No  correction  was  included  in  the  calculations 
for  target  thickness  since  the  thickest  target  was  only  8.7%  of  the 
narrowest  resonance  width. 

Figures  B.16-B.18  show  the  experimental  differential  cross  sections 
for  each  of  the  three  samples.  The  experimental  uncertainities  in 
each  of  the  experimental  points  is  due  to  counting  statistics  and 
the  energy  stability  of  the  machine.  Counting  statistics  lead  to 
errors  as  large  as  ±50%  at  low  values  of  the  cross  section  and  to 
•2%  at  large  values  of  the  cross  section.  The  uncertainities  due 
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Sections 


Due  to  the  large  uncertainities  in  the  experimentally  measured 
cross  section  it  is  impossible  to  determine  which  of  the  theoretical 
curves.  Figures  B.8-B.15,  best  fits  the  experimental  data.  However, 
it  is  possible  to  estimate  the  resonance  energies  and  resonance 
widths  from  the  experimental  data  and  compare  them  with  the  values 
given  by  Aj zenberg-Selove  (164).  These  results  are  tabulated  in 
Table  B.2. 
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TABLE  B.2 

Comparison  of  Experimental  Resonance  Parameters 


Target 

E (keV) 

°1 

riabl(keV) 

E (keV) 
°2 

riab2(keV) 

0.33  keV 

1004  ± 

3 

136  ± 

15 

1211 

± 

1 

21  ± 2 

1.30  keV 

1010  ± 

5 

156  ± 

10 

1213 

± 

1 

24  + 1 

1.95  keV 

999  ± 

15 

161  ± 

7 

1210 

± 

1 

26  ± 1 

Literature 

1028  ± 

10 

140  ± 

10 

1210 

± 

3 

22.5  + 1 

Value  (164) 


•?; 
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